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SECTION  I 
INTRODUCTION 


In  deriving  far  field  conditions  for  unsteady  flows  at  high 
subsonic  Mach  numbers,  one  usually  makes  the  assumption  that  in  the 
distant  field  the  deviations  from  a  parallel  flow  with  the  given 
free  stream  Mach  number  are  small.  Two  quantities  are  of  importance, 
the  ratio  of  the  velocity  perturbations  to  the  free  stream  velocity 
of  sound  and  the  ratio  of  the  velocity  perturbations  to  the 
difference  between  the  free  stream  velocity  and  the  sonic  velocity 
(the  velocity  for  which  velocity  of  sound  and  velocity  of  motion 
are  the  same) .  Accordingly,  it  is  not  permissible  to  choose  a 
boundary  of  the  computed  part  of  the  flow  field  for  which  the  local 
Mach  number  reaches  the  value  of  one.  The  distance  from  the  profile 
for  which  this  prerequisite  is  satisfied  increases  to  infinity  as 
the  free  stream  Mach  number  approaches  one. 

The  original  goal  of  the  research  effort  reported  here  was 
the  derivation,  by  means  of  asymptotic  techniques,  of  far  field 
conditions  for  small  unsteady  perturbations  with  harmonic  time 
dependence  superimposed  to  a  steady  flow  with  the  free  stream  Mach 
number  one.  It  was  planned  to  illustrate  the  mathematical  technique 
by  applying  it  to  problems  with  a  high  subsonic  free  stream  Mach 
number.  This  attempt  led  to  some  developments  beyond  those  found 
in  the  literature.  Usually  one  assumes  that  the  perturbations 
arriving  at  the  outer  edge  of  the  computed  flow  field  can  be 
approximated  by  cylindrical  or  spherical  waves  whose  origin  is 
known.  (Actually,  it  must  be  estimated.)  The  analysis  leads  to 
equations  for  which  the  wave  fronts  may  have  a  more  general 
character,  although  one  must  assume  that  they  are  rather  smooth. 

The  information  provided  in  the  standard  methods  by  estimating 
the  location  of  the  origin  is  then  derived  from  flow  field  data. 
Furthermore,  it  is  a  rather  simple  step  to  proceed  from  perturba¬ 
tions  which  are  harmonic  in  time  to  general  unsteady  perturbations. 

Of  course,  asymptotic  results  are  not  always  applicable.  The 
requirement  that  the  frequency  of  harmonic  oscillations  is  fairly 
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high  or  that  changes  of  the  flow  field  with  time  are  rather 
quick  is  not  always  satisfied.  Moreover,  the  description  of 
the  perturbation  field  by  means  of  smooth  wavefronts  is  not  always 
possible.  Then  one  must  apply  far  field  conditions  (of  a  more 
complicated  nature)  in  which  such  assumptions  are  not  made  (Ref. 

1  and  2) .  Of  theoretical  interest  is  the  relation  of  surfaces 
of  constant  phase  (which  play  an  important  role  in  the  asymptotic 
considerations)  to  characteristic  surfaces  in  the  space  spanned 
by  the  space  coordinates  and  the  time. 

Accordingly,  this  report  consists  of  two  nearly  independent 
parts,  one  consisting  of  Sections  II  to  VI  which  deals  with  far 
field  conditions  at  a  high  subsonic  Mach  number  and  the  remainder 
which  treats  far  field  conditions  for  oscillatory  flows  at  a  free 
stream  Mach  number  one. 

Regarding  the  latter  we  make  the  following  observations. 

We  deal  with  small  unsteady  perturbations  with  harmonic  time 
dependence  superimposed  to  a  steady  flow  with  the  free  stream 
Mach  number  one.  Because  of  the  so-called  freezing  of  the  Mach 
number  distribution  in  the  flow  field  for  free  stream  Mach  numbers 
in  the  vicinity  of  one,  the  results  are  also  applicable  for  flow 
fields  with  such  Mach  numbers.  In  this  study  the  linearization 
is  carried  out  for  the  vicinity  of  a  steady  flow  perturbed  by  the 
presence  of  a  body,  (while  in  the  derivation  for  subsonic  Mach 
numbers,  one  linearizes  for  the  vicinity  of  a  parallel  undisturbed 
flow) .  It  is  assumed  that  the  boundary  at  which  the  far  field 
conditions  are  to  be  applied  lies  at  such  a  distance  that  the  flow 
field  can  be  approximated  by  the  dominant  term  in  the  development 
with  respect  to  distance.  For  these  perturbations  analytic 
expressions  are  available,  even  for  the  axisymmetric  problem.  They 
have  been  found,  simultaneously  by  a  number  of  authors  (Randall, 
Mueller  and  Matschat,  and  Euvrard,  Refs.  3,  4,  and  5).  We  shall 
use  the  form  given  by  Randall.  The  conventional  treatment  of 
far  field  conditions  for  subsonic  flows  makes  use  of  analytical 
expressions  for  the  unsteady  linearized  far  field.  For  flow  fields 
with  a  free  stream  Mach  number  one,  such  expressions  are  not 


available.  One  can  dispense  with,  information  of  this  kind  by 
asymptotic  considerations.  This  is  a  technique  used  in  physics 
to  make  the  transition  from  wave  to  ray  optics.  The  asymptotic 
expressions  are  valid  if  the  frequencies  are  sufficiently  high. 

Such  techniques  have  been  applied  before  for  problems  of  sound 
|  radiation  (for  recent  applications  see  Bayliss,  Gunzburger ,  and 

f  #  Turkel,  Ref.  6,  71.  The  mathematical  nature  of  conventional  f._r 
I  field  conditions  seems  to  indicate  that  one  comes  into  the  realm  of 

Isuch  high  frequency  approximation  as  the  free  stream  Mach  number 
approaches  one . 


# 
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SECTION  II 

THE  SUBSONIC  PROBLEM,  BASIC  EQUATIONS,  PARTICULAR  SOLUTIONS 


Let  x,y  be  a  system  of  Cartesian  coordinates,  t  the  time, 
and  J  the  velocity  potential  which  describes  the  deviations  from 
a  parallel  subsonic  flow.  The  linearized  differential  equation 
which  we  take  as  point  of  departure  for  the  present  investigations 
can  be  obtained  either  by  linearizing  the  potential  equation  for 
unsteady  flow  for  the  vicinity  of  a  steady  parallel  flow,  or  by  a 
transformation  of  the  equation  for  the  propagation  of  sound  in  air 
at  rest.  The  differential  equation  for  the  latter  problem  is 
given  by 


4> _ +  <J>__  -  a"*2<J> _ =  0  (1) 

xx  yy  tt 

where  "a"  is  the  constant  sound  velocity.  The  equation  for 
perturbation  in  air  that  moves  from  right  to  left  with  the 
velocity  U  is  obtained  by  the  following  transformation 

?(£,y,t)  =  $(*,?,£) 


where 

x  =  x  +  Ut 

y  -  y  (2) 

t  ■  t. 


One  obtains 

(l-U2/a2)$~~  +  -  2 (U/a2) -  (l/a2>$£€  =0.  (3) 
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We  introduce  the  free  stream  Mach  number 


M  ■  U/a 

make  all  length  dimensionless  with  a  characteristic  length  L,  for 
instance  the  chord  of  the  profile, 

x  =  x/L  ,y  =  y/L 

and  introduce  a  dimensionless  time 

t  =  a  E/L  =  a  t/L.  (4) 

Then  one  has  instead  of  Eq.  (2) 

x  =  L  (x-M  t). 

With  $<x,y,t)  =<jj(x,y,t),  one  then  obtains  from 

(1"M2)^xx  +  *yy  “  2M*xt  "  *tt  * 

This  derivation  emphasizes  the  fact  (actually  a  rather  obvious 
one) ,  that  perturbations  once  introduced  behave  like  perturbations 
in  still  air,  except  that  they  float  downstream  with  the  velocity 
U.  In  the  last  equation,  U  is  replaced  by  the  dimensionless 
velocity  M.  We  now  consider  particular  solutions  which  are 
periodic  in  time  with  the  circular  frequency  v 

$(x,y,t)  ■  <Mx,y)  exp(ivE) 

or 

$<x,y,€)  *  <j> (x,y)  exp(i (vL/U)  (Ut/a) ) 

*  <J>(x,y)  exp(iy'Mt)  (6) 

*  <j>(x,y)  exp(iwt) 


Eq.  (3) 

0  (5) 
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Mil a** 


In  the  second  version  of  Eq.  (6)  the  familiar  reduced  frequency, 
here  denoted  by  p ' ,  has  been  introduced 

U*  -  (vL/U)  (7) 

Furthermore,  we  have  set 

a  =  p'M  =  (vL/a)  (8) 

One  then  obtains 

( 1-M2 )  <pxx  +  <j>  -  2iuM<(,x  +  aj2(j>  =0  (9) 

In  Ref.  5  Eq.  (9)  has  been  transformed  in  the  following 
manner.  First  one  carries  out  the  Prandtl  Glauert  coordinate 
distortion 


x  =  x 


„2. 1/2 

y  =  (1-M  )  y 


One  obtains 

^xx  +  ^yy  ~  2iu>M(l-M2)  +  (1-M2)  ^u)24>=  0 

Next  one  removes  the  term  with  <}>£  by  the  transformation 
<P(x,y)  =•  exp(iu>M(l-M2) -1x)  $  (x,y) 

Introducing 

U  *  (1-M2)_1U  =  Md-M2)*1!!' 


(10) 


(ID 


(12) 


(13) 


One  obtains 


yy 


2* 
V  y 


0 


(14) 
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This  is  the  Helmholtz  equation.  Usually  the  Helmholtz  equation 
arises  in  the  study  of  periodic  perturbations  in  still  air.  One 
then  substitutes  the  following  expression  into  Eq.  (1) 

J(x,y,T)  •  J (x/L)  ,  (y/L)  )exp(ivt) 


Using  Eq.  (8)  one  obtains 


s27  + 

3  (x/L)  2 


d2j 

3  (y/L)  2 


2 — 

0)  (j) 


0. 


(15) 


Notice  that  Eqs.  (14)  and  (15)  have  different  factors  of 
4>  and  ?.  The  substitute  frequency  u  tends  to  infinity  as  M  tends 
to  one.  This  suggests  that,  in  the  vicinity  of  Mach  number  one, 
the  high  frequency  limits  of  the  far  field  conditions  is  applicable. 
The  transformations,  Eqs.  (10)  and  (12)  ,  amount  to  a  rather  com¬ 
plicated  distortion  of  the  original  flow  pattern.  There  is  no 
direct  physical  relation  between  the  flow  fields  described  by  the 

A  - 

functions  <t>  and  <p . 

Particular  solutions  of  Eq.  (14)  are  given  by 


Here 


♦m(x,y)  = 


H(2)  (ur-\  fCos(me) 
Hm  {Ur)tsin(m6) 


r  =  [Si2  +  $2]1/2  =  [x2  +  (1-M2)y2)1/2 


(16) 


(17) 


COS  0 
sin9 


x/r  =  x/r 

y/r  =»  (l-M2)1/2y/r 


0  <  e  <  2tt 


(18) 


(2) 

Hm  is  the  Hankel  function  of  order  ra,  where  m  is  a  positive 
integer  or  zero. 

One  has  for  large  values  of  a  complex  variable  z 
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const  z 


(19) 


) exp(-iz) 


where  P  stands  for  a  semiconverging  power  series  in  z**1.  The  exact 
form  of  the  constant  is  unessential  in  the  present  context.  These 
particular  solutions,  partially  written  in  terms  of  the  original 
coordinates ,  read 


<p  -  exp(i<uM(l-M2) -1x) 


H^UU-mV^x2  .  <1-M2)y2)1/2) 


(20) 


and  asymptotically 

<p  =  const  r~1/2exp  {iw (1-M2)  _1  [Mx  -  (x2+(l-M2)y2)  1/2\ 

2  —1 

The  factor  (1-M  )  in  the  exponential  function  and  in  the  argument 
of  the  Hankel  functions  makes  it  obvious  that  these  solutions  will 
fail  as  the  Mach  number  approaches  one. 

An  intuitive  picture  can  be  obtained  in  the  following  manner. 
Assume  that  one  generates  perturbations  by  short  pulses  spaced  at 
equal  time  intervals  in  a  flow  with  Mach  number  M.  We  consider  the 
wave  front  pertaining  to  each  perturbation  at  some  later  time  t. 

We  mentioned  previously  that  these  fronts  spread  out  as  in  air  at 
rest  (that  is,  with  the  sound  velocity  "a")  but  at  the  same  time 
move  downstream  with  the  dimensionless  velocity  M.  For  a  Mach 
number  smaller  than  1,  the  system  of  such  wave  fronts  is  drawn  in 
Figure  1.  The  distance  between  two  adjacent  wave  fronts  is  smaller 
on  the  upstream  side  of  the  point  where  the  perturbations  are 
introduced  and  larger  at  the  downstream  side  compared  with  the 
wave  fronts  in  air  at  rest  (which  would  be  formed  by  concentric 
circles) .  At  a  free  stream  Mach  number  1  the  upstream  distance 
reduces  to  zero  (Fig.  2).  This  accumulation  of  perturbations 
suggests  the  difficulties  which  arises  in  a  flow  with  the  free 
stream  Mach  number  1.  It  is  obvious  that  in  such  a  flow  the  local 
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(21) 


deviations  of  velocity  of  sound  and  of  the  particle  motion  from 
those  of  the  free  stream  conditions  can  no  longer  he  disregarded. 
A  discussion  where  this  modification  is  made  will  he  carried  out 
later. 
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SECTION  III 


THE  SUBSONIC  PROBLEM,  HIGH  FREQUENCY  APPROXIMATION 

To  obtain  a  high  frequency  approximation  to  solutions  of 
Eq.  (9)  we  make  the  transformation 

<j>(x,y)  =  h(x,y)  exp(-iujg(x,y) )  (22) 

and  satisfy  Eq.  (9)  in  the  dominant  powers  of  w. 

One  has 

♦x"I(Vh)  “iuj9x]<i> 

$Xx  =  [  (hxx/h) -2  (hx/h)  itogx  -  u>2g2  -  iwg^]  <j> 

<J>yy  =  1  ^yy^  “2  (hy/h)  ift,gy  “  **  ^yy^* 

Then  from  Eq.  (9) 

-u2 [ (1-M2) g2  +  g2  +  2Mgx  -  1] 

-iuh"1[(l-M2)2hxgx  +  2hygy  +  2Mhx+  h((l-M2)gjot  +  g^)  ]  (23) 

+h_1t(l-M2)hxx  +  h^  ]  =  0. 

2 

The  dominant  terms,  (those  with  the  factor  w  )  determine  the 
function  g(x,y) .  One  obtains  the  equation 

(1-M2) g2  +  g2  +  2Mgx-l  -  0.  (24) 

The  function  h(x,y)  depends  upon  u.  Including  one  term  beyond  the 
lowest  order  approximation  we  set 

h(x,y,w)  -  H1(x,y)+ui'*1H2(x,y)  (25) 
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The  terms  of  order  <u  and  of  order  1  in  Eq.  (23)  give  respectively 


Hl,xIa'H2)gx+MI  +  Hl,y9y  +  U/2)H1[(tM2)gxx+gyy]  -  0  (26) 

H2(J!ta-M2)gx+Ml  +  H2(ygy  +  (1/2)H2I(1-M2)gxx  +  gyyJ 

(27) 

•Hi/2)[(l-M2)HljXX  +  Hl  yy]  -  0. 

Eq.  (24)  is  a  first  order  partial  differential  equation  in  one 
dependent  variable.  Applying  a  standard  procedure  one  differenti¬ 
ates  Eq. (24)  with  respect  to  x  and  y 

[U-M2)gx  +  Ml  9  xx  +  Vxy  =  ° 


and 

Ul-M2)gx  +  M]gxy  +  gygyy  =  0. 


We  denote  by  D  differentiation  in  a  direction  given  by 


(dy/dx)  =  gy/[(l-M2)gx  +  M] 


Let  Dl  be  the  line  element.  Then  one  has 


Dx 
DT  = 


(X-M2)gx  +  M 
(g2  +  ((1-M2)gx+M)2]1/2 

_ £y _ 

Igy  +  ( ( l-M2 )  gx+M)  ^  ] 1/1 2 


D 

Dl 


DX  _3 

51  3x 


+ 


Dy  3 
"51  Sy 


(28) 


(29) 


(30) 


(31) 
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Curves  with  the  slope  given  by  Eq.  (29)  are  called  characteristics. 
Then  from  Eqs.  (28) 


Dg^Di  =  0 

and  (32) 

Dgy/Di  =  0 . 

Furthermore 

Dq/Dl  =  gx(Dx/Di)  +  gy(Dy/D£)  (33) 

and,  from  Eqs.  (26)  and  (27) 

DdogH^/Di  +  (1/2)  [  (1-M2)  g^  +  9^1/lg*  +  ((l-M2)gx+M)2]1/2=o 

DH2/D£+(l/2)H2[(l-M2)gxx+gyy]/[g2  +<(1-M2)gx  +  M)2]1/2  (34) 

+  (i/2) I(l-M2)Hlxx+H1>yy]/[g2  +  ((1-M2)gx  +  M)V/2  -0 

It  follows  from  Eqs.  (32),  that  along  the  characteristics  . 

g  ■  const  and  g  ■  const  and  subsequently,  from  Eq.  (29) ,  that 
x  y 

these  curves  are  straight  lines. 

For  a  preliminary  discussion,  we  consider  characteristics 
which  st:art  at  the  origin.  Then 

Dy/Dx  *  y/x 


and  because  of  Eq.  (29) 

gy  ■  (y/x)[(l-M2)  gx  +  m] 

One  obtains,  by  substituting  this  expression  into  Eq.  (24) 
and  by  solving  the  resulting  quadratic  equation. 
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(35) 


9x  =  +  x/(x2  +  (1-M2)y2)  1/2J 

gy  *  y/(x2  +  (1-M2)y2)1/2 
Moreover/  if  g  =  0  for  x  *  0  and  y  *  0 


g  =  gxx  +  gyy 

g  -  (1-M2)_1[-Mx  +  (x2  +  (1-M2)y2)1/2]  (36) 

Consider  acurveg(x,y)  =  K  (where  K  is  some  constant).  One  finds 
from  the  last  equation  that  such  a  curve  satisfies 

(x  -K  M)2  +  y2  -  K2.  (37) 

This  is  a  circle  with  a  radius  |k|  and  center  at  the  point 
(x,y)  =  (KM,0) .  For  K  positive  these  are  the  circles  shown  in 
Figures  1  and  2. 

The  limit  M  *  1  is  readily  discussed.  One  obtains 


gy  *  Y/ |*| 

gx  =  (I-mWm  +  (x/|x|)  ]  . 


(38) 


Hence 


gx  =  (1+M) 


-1 


1/2,  x  >  0 


lim  _  lim  -1  ^ 

M-l  9x  *  M-l  "{1_M)  “  x  <  0* 


Actually,  for  M  ■  1  no  perturbations  can  arrive  at  stations  x 


(39) 
<  0. 
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Eq.  (37)  shows  that  the  factor  exp  (-iwg(x,y))  in  Eq.  (22) 
anticipates  the  waviness  of  the  flow  field.  The  function  h(x,y) , 
gives  a  modulation  of  the  amplitude  and  since  it  may  be  complex 
also  some  modification  of  phase.  The  phase  modulation  is  expected 
to  be  small.  The  evaluation  of  is  best  carried  out  on  the  basis 
of  Eq.  (26)  rather  than  Eq.  (34) .  One  obtains  for  the  coefficient 
of 

(1/2)((1-M2)gxx  +  gyyI  =  (1/2)  r"1 


where  according  to  the  definition  in  Eq.  (17) 


r  *  [x2  +  (1-M2)y2]1/2, 


Using  the  second  of  Eqs.  (35)  and  Eq.  (29)  ,  one  obtains  from 
Eq.  (26) ,  after  multiplication  by  r  H^-1 

(Hl,x/Hl)x  +  (Hl,y/Hl)y  +  (1/2)  =  °* 

This  equation  is  solved  by 


H  *  f(x/y)  [x2+y2r1/4.  (40) 

Using  Eq.  (22)  and  Eq.  (36)  one  thus  obtains  the  following 
approximation  for  <j> 

<J>  *  f  (x/y)  (x2+y2)  -1//4exp(-iwg(x,y) 

(41) 

=>  f  (x/y)  (x2+y2)"1/4exp[iu>(l-M2)_1(Mx  -  (x2+(l-M2)y2) 1/2]  . 

This  is  the  first  term  in  the  asymptotic  expressions  for  exact 
solutions  shown  in  Eq.  (21).  Comparing  Eqs.  (41)  and  (21)  one 
obtains 


f(x/y)-l(x2  +  y2)/(x2+(l-M2)y2]1/4 


cos  (mO) 
sin (m6 ) 


(42) 


where  6  is  defined  in  Eq.  (18) 


F 


The  asymptotic  expression,  Eq.  (21) ,  arises  from  the 
development  of  the  Hankel  functions.  The  argument  of  the  Hankel 
functions  tends  to  infinity  even  as  M  tends  to  1  for  finite  x  and 
y.  This  justifies  the  high  frequency  approximation.  The  defini¬ 
tion  for  8  in  Eq.  (18)  shows  that  for  M-*l,  y  finite  and  positive 
and  x  >0,  e  tends  to  zero,  and  for  M-*-l,  y  finite  and  x  <  0,  8 

tends  to  tt.  For  M  close  to  1  the  function  f(x/y)  therefore  changes 
very  rapidly  in  the  vicinity  of  x  =  0 .  This  corresponds  to  the 
accumulation  of  perturbations  shown  in  Fig.  2.  Notice  also  that 
the  function  f(x/y)  has  a  strong  peak  at  x  =  o  for  M-+-1  and  assumes 
the  value  (1-M2)  “1/4 . 

In  preparation  for  a  generalization,  we  rewrite  these  equations 
in  terms  of  the  direction  normal  to  the  wave  fronts  (given  by 
curves  g  *  const)  and  in  terms  of  the  radius  of  curvature.  It 
was  shown  above,  that  lines  g  =  const  are  circles  which  float 
downstream  with  the  dimensionless  velocity  M;  at  time  t  the  circles 
have  the  dimensionaless  radius  R  =  t.  Let  a  be  the  angle  of  the 
normal  to  the  wave  front  with  the  x-axis  at  some  point  of  the  wave 
front  (Fig.  3)  and  let  8  be  the  angle  of  the  ray  AB  with  the 
x-axis.  Then  one  has 


x  =  R(M  +  cosa) 
y  =  R  sina 

R1  =  (x2  +  y2)1/2  -  R(M2+2Mcosa  +1) 1/2 

tg  8  *  y/x  =  sin  a/(M+cosa) 

2  -1/2 
sinB  =  sina  (M  +2M  cosa  +1)  ' 

cos 8  =(M  +  cosa)  (M2  +  2M  cosa  +  1)"1//2 


(43) 


In  Eqs.  (17)  and  (18)  the  quantities  r  and  6  have  been  defined. 
One  obtains 


r  -  R(l+M  cosa) 
tg  8  -  (i-M2)1/2tg  8 


(44) 


The  expression  g,  Eq.  (36) ,  reduces  to 

g  =  R 

With  a  modified  function  f,  one  then  obtains  from  Eq.  (41) 

4  «  f  (a)R“1/2exp(-iu,R)  (45) 

The  direction  of  the  characteristic  is  given  by  a  line  Q  =  const, 
or 

Dx/Dy  =  sina/(M  +  cosa)  (46) 

To  visualize  what  happens  in  a  more  general  case,  in  which 
the  initial  line  g  =  const  is  not  a  circle  but  a  general  closed 
curve,  one  uses  the  idea  that  lines  g  =  const  at  a  later  time 
arise  from  an  initial  line  in  the  same  manner  as  in  air  at  rest, 
except  that  they  float  downstream  with  the  dimensionless 
velocity  M.  We  already  found  that  the  characteristics  are  straight 
lines  in  any  case.  For  a  point  on  a  line  g  =  const  with  normal 
given  by  a>  the  slope  of  the  characteristic  is  then  given  by 
Eq.  (46). 

The  above  description  suggests  that,  along  a  specified 
characteristic,  <j>  is  determined  by  the  curvature  of  the 
initial  line  g  =  const  in  the  same  manner  as  for  perturbations 
emanating  from  a  circle.  This  motivates  the  following  representa¬ 
tion;  let  the  inital  curve  g  =  const  be  given  by 

x  = 

y  = 

where  s  is  the  arc  length 

(dxQ/ds) 2  +  (dyQ/ds) 2  *  1.  (48) 
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Then  the  angle  of  the  normal  to  this  curve,  denoted  again  by  a, 
is  given  by 


cosa  =  dyQ/ds 
sina  =  -(dxQ/ds) 


one  has 


gx  =  jgrad  gjcosa 
9y  ■  | grad  g|sina. 

By  substituting  the  above  equations  into  Eq.  (24)  one  obtains 

(grad  gj  =  1/(1+M  cosa) 


(49) 


and  therefore 


g  *  cosa  /(1+M  cosa) 

(50) 

g^  *  sina  /(1+M  cosa) 

The  direction  of  the  characteristics  is  given  by  Eq.  (29) .  One  has 
for  the  denominator 

(1-M2)gx  +  M  *  (M+  cosa) /(1+M  cosa) 

Therefore 


Dy/Dx  -  sina/(M+  cosa)  (51) 

With  a  parameter  p  one  then  obtains  the  following  parametric 
representation  for  the  characteristic 
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(52) 


x  =  x  +  (p-p  )  (M  +  cosa) 
o  o 

y  =  yo  +  (p-pQ)sina 

We  have  denoted  the  arc  length  of  a  characteristic  by  d l.  One 
obtains  from  the  last  equation 

dil/dp  =  (M2  +  2M  cosa  +  1)1/2.  (53) 

Furthermore , 

dg/dp  =  g  (dx/dp)  +  g  (dy/dp)  =1.  (54) 

x.  y 

It  is  shown  in  Appendix  II,  that  p  can  be  identified  with 
the  radius  of  curvature  of  the  line  g  =  const,  to  be  denoted  by  R. 
The  function  H1  is  determined  from  Eq.  (26)  .  Using  Eq.  (50)  and 
Eq.  (34)  one  obtains 

(55) 

(d(logH,)/dJl)  +  (1/2)  [  (1-M2)  g  +g  ]  ( 1+M  cosa)  (1+2M  cosa  +  M2) -1/2=0 
±  xx  yy 

In  Appendix  II  it  is  shown  that  this  leads  to 

d(log  H^/dp  +  (1/2)  (1/p)  =  0  (56) 

Hence 

H1  =  f x (s)  p“1/2 

For  circular  waves  p  is  given  by 

p  -  R  ■  Rj/fl  +  2M  cosa  +  M2)1//2  (57) 

in  other  words,  p  is  proportional  to  R1  with  a  factor  of 
proportionality  which  is  different  for  different  characteristics. 

It  is  also  shown  that 


H2  -  f2(s)p"1/2. 

In  the  derivation  of  the  expression  for  H2  in  Appendix  I I 
the  assumption  has  been  made  that  p  is  large  (while  p^  is  not  too 
large) .  The  last  formulae  therefore  gives  a  good  approximation 
only  at  a  sufficient  distance  from  the  initial  curve.  The 
parameter  s  characterizes  the  individual  characteristics . 

We  shall  use  Eq.  (22)  in  the  asymptotic  form  of  Eq.  (41) , 
but  somewhat  extended 

♦  =  [f1  (s)  Rj^172  +  f2(s)R1“3/2]exp(-iwg)  (58) 

With  Eq.  (57)  and  modified  functions  f^  and  f2  one  can  also  write 

<J>  =  (f1(s)  p-1^2  +  f2(s)p”3//2]exp(-iu>g)  .  (58') 


SECTION  IV 


ASYMPTOTIC  FAR  FIELD  CONDITIONS 


Asymptotic  far  field  conditions  are  obtained  from  Eq.  (58')  in 
the  following  manner.  One  forms  the  derivative  in  the  direction 
of  the  characteristic.  Then  no  derivatives  of  the  unknown 
functions  f^(s)  and  f^Cs)  will  appear.  Using  Eq.  (54)  one  obtains 

D4>/DH  =(Dp/DJL)  [f x (s)  (-(l/2)p“3/2-iwp“1/2)+f2(s)  (-(3/2)p"5/2-iop"3/2)  )exp(-iwg)  (59) 


The  expression  $  of  Eq.  (58')  then  satisfies 

lir  +  if  <^  +  2|(5i+  i„)i*  =  0.  (60) 

The  second  operator  (which  is  applied  first)  makes  the  term 
f1(s)p  1//,2(-iajp)  zero  and  changes  the  expression  f2  (s)  p”3//2exp(iajg) 
into  -  dp/dJl  f2  (s)  ~5//2exp(-iug)  .  The  first  part  of  the  operator 
(which  is  applied  afterwards)  causes  this  term  to  vanish.  For 
large  values  of  p  (or  of  u)  one  may  disregard  the  term  p  '  and 
then  one  obtains  the  simple  far  field  condition 


(J, 

DJI  '2p 


+  iw)  ]  <j>  *  0 


(«•) 


Also,  the  term  l/(2p)  may  be  omitted  if  u>  is  sufficiently'  large 


(62) 


The  last  formulation  corresponds  to  the  original  formulation  of 
Sommerfeld.  Eqs.  (61)  and  (62)  are  different  forms  of  the  far 
field  conditions  of  Bayliss,  Gunzberger,  and  Turkel  (Ref.  6,7). 

The  expression  Eq.  (60)  written  in  detail  gives 


(D£)2 

lDjT 


+ 


3 

P 


iw-u)  ]  <)> 


=0. 


(63) 
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The  occurrence  of  a  second  derivative  in  the  boundary  conditions 
for  a  second  order  differential  equation  is  unusual.  In  the 
present  problem  second  derivatives  cannot  be  avoided ,  but  they 
can  be  transformed  by  means  of  the  differential  equation  into  second 
derivatives  along  the  contour  of  the  computed  region.  If  one  uses 
higher  order  boundary  conditions  of  the  Bayliss,  Gunzberger,  Turkel 
type,  then  even  higher  derivatives,  first  in  the  direction  of  the 
characteristics  and,  after  a  transformation,  along  the  contour  of  the 
computed  region  are  encountered .  Exact  boundary  conditions  for 
the  problem  have  been  derived  in  Ref.  5;  they  are  of  a  global  nature 
(that  is,  all  boundary  points  interact  with  each  other) .  The 
occurrence  of  higher  derivatives  means  that  the  solution  in  some 
vicinity  of  the  point  under  consideration  enters  the  formulation. 

In  this  manner  the  underlying  global  character  of  such  far  field 
conditions  makes  itself  felt. 

Some  further  details  are  added.  Let  3  be  the  angle  of  the 
direction  of  the  characteristic  with  the  x  axis.  Then,  according 
to  Eq.  (43) , 


sin  0  =  sin  a  (M2+  2M  cos  a+  1)  1//2 

2  -1/2 
cos  0  ■  (M  +  cosa) (M  +  2M  cos  a+  1)  ' 


(64) 


(By  our  definition  a  is  the  angle  of  the  normal  to  the  line  g  * 
const  with  the  x-axis.)  It  may  be  desirable  to  express  various 
quantities  in  the  terms  of  the  angle  6 .  We  introduce  an  auxiliary 
anglefi  (Fig.  3). 

Then 

sin  5s  M  sin  0 

a  ■  (0+5)  (65) 

(M2  +  2M  cos  a+  1)  ^2  ■  M  cos  0  +  cos  6 
DR/D&  ■  Dp/Di  *  (M  COS  0+  COS  6)”1 


t 
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2  2 

Now  we  express  the  derivative  d  $/d l  .  We  introduce  a  local 
system  of  coordinates  i,  where  C  has  the  direction  of  the 
characteristic,  n  the  direction  of  the  tangent  to  the  boundary 
of  the  computed  region.  Then  d/di  =  d/d£.  Let,  at  the  point 
under  consideration,  y  be  the  angle  of  the  n  axis  with  the  x 
axis  (Figure  4) . 

Then 

X  -  XQ  ~  £  COS  $  +  T)  COS  Y 

y  -  yQ  =  K  sin  3  +  n.  sin  y 

and 

2  2 
=  ♦xx  COS  6  +  *xy  Sin  23  +  ^yy  Sln  B 

4)nn  =  <t>xx  cos2Y  +  ♦  Sin  2y  +  (fyy  sin2Y  (66) 

"  ^xx  cos  8  cos  Y  +  4>xy  sin(S+Y) +<J>yysin  3  sin  y. 

In  addition,  one  has  the  differential  equation  for  <J>  (Eq.(9)) 

2io>M<j>x  -  U)2<j>  =  (1-M2)<j>xx  +  ^yy * 

Eliminating  4*^,  ^yy'  and  4>xy  ^rom  these  four  equations  one 
obtains 

<(>££  [l-M2sin2Y]+4>nr)  (l-M2sin2B)  [2cos  (B-y)  +M2  (cos (B+y)  -cos (B~y)  )  1 

(67) 

-(2iwM<j>x  -  w  <j>)  sin^  (B-Y)  -  0. 

This  equation  is  verified  in  Appendix  III.  It  allows  one  to 
2  2 

express  d  4>/diL  =  in  terms  of  derivatives  of  4>  and  <p^  xn  the 
n  direction  (that  is  in  the  direction  of  the  contour  of  the  computed 
part  of  the  flow  field) . 


In  practice,  one  will  choose  the  boundary  of  the  computed 
part  of  the  flow  field  so  that  the  n  direction  coincides  with 
either  the  x  or  the  y  direction.  Then  considerable  simplification 
will  occur. 

These  results  can  be  applied  in  two  different  forms.  In  a 
procedure,  which  corresponds  to  that  of  Bayliss,  et  al.,  one 
assumes  that  all  waves  originate  at  the  same  origin  (somewhere  in 
the  middle  of  the  profile) .  The  direction  of  the  characteristics 
is  then  given  by  rays  through  the  origin.  The  simplest  form  is 
obtained  from  Eq.  (58) 

<j>  =  [^(s)^"172  +  f2(s)R1“3/2]exp('-iu)g)  . 

where 


Rx  =  (x2  +  y2) 1/2  (68) 

and  according  to  Eq.  (36) 

g  =  (l-M^^t-Mx  +  (x2  +  (1-M2)y2)1/2)  .  (69) 

According  to  the  derivation  of  Eq.  (36)  one  has 

Dg/Di-  gx  (Dx/DJl)  +  gy  (Dy/DJl)  =  gx(x/R1)  +  g^y/R^ 

Dg/D£=  g/Rr 

According  to  Eq.  (54)  ,  one  has  dg/dp  *=  1,  Therefore,  from  Eq.  (62) 

(D0/DJI)  +  iwtg/R^fj)  =  0  (70) 


and  from  Eq.  (61) 

(D(j)/D£)  +  [iwfg/R^  +  (1/(2R1))J^  *  0 
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The  counterpart  to  Eq.  (63)  is 

(D24>/D*2)  +  ( ( 3/R^  +  ( 2iwg/R1)  ]  (D<f>/DA)  +  [  ( 3/4R2)  +  ( 3iu>g/R2)  -  (gw/R^ 


(71) 

2R=0 


A  second  form  of  the  far  field  condition  arises,  if  one 
determines  the  likely  origin  for  the  waves  which  arrive  at  some 
point  of  the  contour  of  the  computed  flow  field  from  the  function 
<t>  and  its  derivatives  along  the  contour.  One  proceeds  in  the 
following  manner.  The  analysis  is  based  on  the  assumption  that 
the  function  h  changes  only  slowly.  In  principle,  h  and  <f>  are 
allowed  to  be  complex.  Then  one  has 

arg(<j>)  =  arg(h)  -  aig 


and 


g  = 


w  1[arg(h) 


arg(<j>)  ] 


locally  arg(h)  can  be  replaced  by  a  constant. 

From  this  expression  one  obtains  gx  and  gy  along  the  contour 
of  the  computed  region  by  numerical  differentiation,  and  hence 


cos  a  -  gx/(g2  +  9y)1/2 

(72) 

sin  a  =  gy/(g2  +  9p1/2 

As  a  check  (which  may  lead  to  small  corrections)  one  has  from  the 
first  of  Eqs.  (50) , 


|  grad  g|  =  (1  +  M  cosa) 


The  fact  that  g  =  p  is  identical  with  the  radius  of  curvature  of 
a  line  g  =  const,  is  hard  to  apply  in  a  practical  computation, 
because  one  needs  interpolations  to  identify  these  lines. 


In  Appendix  II  the  following  formulae  have  been  derived 

-1  -3  2 

g^  =  R  {1  +  Mcos  a)  sin  a 

g^,  =  R*~^(l  +  M  cos  a)  2  (M  +  cos  a)  2 

g  *  R  1(1  +  Mcos  a)'3(-sin  ot(M  +  cos  a)), 
xy 

These  equations  can  be  rewritten  as 

-1  -2  2 

gxx/|grad  g |  -  R  (1  +  M  cos  a)  sin  a 

gyy/|grad  g |  *  R-1(l  +  M  cos  a)”2(M  +  cos  a)  2 

-1  -2 

9xy/,|grad  gl  "  R  +  M  cos  a)  (-sin  a(M  +  cos  a))  • 

Assume  that  the  boundary  of  the  computed  region  is  given  by  a 
line  x  *  const.  Along  this  line  one  can  determine  g  and  g 
(if  one  allows  the  use  of  first  derivative  in  the  direction  of  the 
normal  to  this  line.  Thus,  one  obtains  two  expressions  for  R 
(which  in  some  way  must  be  reconciled  with  each  other) .  If  one 
uses  only  derivatives  along  the  boundary,  then  one  has  only  gyy 
at  one's  disposal. 

Equivalent  to  this  formulation  is  the  following  approach. 

One  determines  along  the  boundary  line  the  value  of  a  from 
Eq.  (72) .  Assume  again  that  the  boundary  is  a  line  x  =  const. 

One  can  then  form  3a/ 3y  numerically .  (Since  a  has  a  direct 
geometrical  meaning,  one  will  be  able  to  judge  whether  a  is 
sufficiently  accurate  and  smooth.)  The  radius  of  curvature  is 
then  obtained  in  the  following  manner.  If,  as  in  Eqs.  (47)  etc., 
ds  is  the  line  element  along  the  contour  g  =  const,  then  one  has 

R-1  =  da/ds  =  -  a  sin  a  +  a  cos  a 

x  y 
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M 


Furthermore  since  a  =  const  along  a  characteristic 


Hence 


a  (M  +  cos  a)  +  a  sin  a  =  0 . 
x  y 


R  ^  =  ay(l  +  M  cos  a)  (M  +  cos  a)  ^ 


R  1  =  “axd  +  M  cos  a)  (sin  a)-1 


In  these  formulae  (R  A)  is  obtained  by  numerical  differentia¬ 
tion  along  the  contour  of  the  computed  region.  Equations  (73)  and 
(74)  are  basically  the  same;  this  is  shown  in  the  following  manner 

a  =  arctg  (gy/gx) 


“y  = 


g  g  -  g  g 
YY  x  ^xy^\ 

(9 x  +  4  > 


^ cos  a-  g  sm  a 


| grad  g 


=(g  cosa-g  sina) (1+Mcosa) 
yy  xy 


Then  from  Eq.  (74) 

“1  2  —7 

R  ~(9vVCOS  a  -  sin  «)  (1  +  M  cos  a)  (M  +  cos  a) 
yy  *y 

Substituting  here  g  and  g  from  Eq.  (73)  one  obtains  indeed  an 

yy  *y 

identity. 

In  the  last  equation  one  can  eliminate  the  mixed  derivative. 

From  the  fact  that  along  a  characteristic  g  and  g  are  const, 

x  y 

one  obtains 


gxxcos  6  +  gxy  sin  3  *  0 


gxycos  8  +  9yy  sin  6  "  0 


hence,  with  Eq.  (64) 
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r 

g  *  -g  sin  a(M  +  cos  a) _A 
xy  yy 

and 

.  2  -2 

gxx  *  gyy  sin  a(M  +  cos  a) 

We  make  the  following  observation.  The  first  form  of  the 
far  field  conditions,  given  by  Eqs.  (62)  through  (65),  can  be 
derived  from  the  asymptotic  development  of  the  Hankel  functions, 

Eq.  (19) .  In  practice,  one  uses  only  one  or  two  terms  of  the 
development  with  respect  to  z~^ .  The  representation  for  fixed  rw 
deteriorates  with  m  (whose  meaning  can  be  recognized  from  Eq.  (16) . 

The  asymptotic  far  field  conditions  (in  either  form)  therefore  are 
applicable  only  at  the  outer  contour  where  the  amplitude  of  <j> 
does  not  change  too  rapidly  and  if  ru  is  not  too  small. 

The  far  field  condition  in  the  original  form  of  Bayliss,  et.  al , 
require  that  one  first  carry  out  the  transformation  (Eq.  (12));  in 
the  present  formulation  one  uses  directly  the  original  form 
(Eq.  (9))  of  the  partial  differential  equation.  Whether  it  is 
worthwhile  to  carry  out  the  transformation  Eq.  (.12)  ,  also  in  the 
region  close  to  the  profile  where  the  differential  equation  is 
more  complicated,  depends  upon  practical  considerations.  The 
two  formulations  are  equivalent. 

In  one  formulation  of  the  exact  far  field  conditions  (which 
have  global  character)  (Ref.  1)  one  replaces  the  effect  of  periodic 
perturbations  in  the  computed  field  by  unknown  periodic  perturba- 
tions  at  the  contour  of  the  computed  flow  field.  The  local 
intensity  of  these  perturbations  is  one  of  the  unknowns  of  the 
problems.  The  field  at  a  fixed  point  of  the  outer  contour  of 
the  computed  flow  field  appears  as  a  superposition  of 
,  perturbations  at  all  other  points  of  the  contour.  This  is,  of 
course,  a  much  more  general  approach  than  the  asymptotic  form, 
which  requires  that  the  perturbations  behave  as  if  they  came  from 
one  point.  The  potential  at  some  point  of  the  contour  is  then 
expressed  by  an  integral  which  contains  contributions  of  the 
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entire  contour.  The  normal  derivative  is  expressed  by  such  an 
integral  but  in  addition ,  it  contains  a  term  which  depends  upon 
the  local  strength  of  the  perturbations.  In  this  formulation  the 
normal  derivative  (rather  than  the  derivative  in  the  direction  of 
the  characteristics)  appears. 

For  the  purpose  of  determining  far  field  conditions,  the 
idea  of  a  nonreflecting  wall  has  been  put  forth.  If  one  uses  a 
local  formulation  of  the  boundary  conditions,  a  wall  can  be  made 
nonreflecting  only  for  waves  of  a  selected  direction.  The 
condition  of  Bayliss  et  al . ,  gives  a  wall  which  is  nonreflecting 
for  waves  that  come  from  the  origin.  In  the  author's  opinion 
this,  rather  than  the  direction  normal  to  the  wall,  is  the 
appropriate  choice.  It  is  worth  noting  that  in  this  formulation 
the  direction  of  the  wall  enters  only  if  one  uses  approximations 
of  a  higher  order. 


SECTION  V 

GENERAL  UNSTEADY  PERTURBATIONS  IN  A  SUBSONIC  PARALLEL ' FLOW 


Equations  (60)  through  (62)  can  serve  to  derive  asymptotic 
far  field  conditions  for  more  general  unsteady  perturbations  in  a 
subsonic  flow,  if  the  changes  with  time  are  rather  rapid.  This 
generalization  is  desirable  for  the  following  reason. 

Upstream  of  the  profile  where  the  waves  propogate  againsts 
the  oncoming  flow,  they  move  fairly  slowly.  Under  these  circum¬ 
stances  the  dependence  of  the  speed  of  motion  and  of  the  local 
sound  speed  upon  the  amplitudes  of  the  waves  become  important. 
Portions  of  the  waves  where  the  pressure  is  higher  travel  faster 
than  the  average.  This  leads  to  a  distortion.  Even  if  at  the  point 
where  they  are  generated  (say  at  the  profile)  the  waves  are 
sinusoidal  in  time,  they  lose  this  property  as  they  travel  over 
some  distance  (although  they  remain  periodic) .  This  effect  is  not 
taken  into  account  in  a  linearized  approach  where  the  unsteady 
perturbations  are  considered  as  small;  there  the  speed  of  wave 
propagation  is  solely  determined  by  the  properties  of  the  under¬ 
lying  steady  field.  The  wave  amplitudes  for  which  the  linearized 
approach  is  sufficient  become  smaller  as  the  Mach  number  approaches 
one.  This  state  of  affairs  is  clearly  seen,  if  one  studies  the 
propagation  of  one-dimensional  waves  in  a  tube.  Accordingly,  it 
may  be  desirable  to  take  nonlinear  effects  into  account  within  the 
computed  part  of  the  flow  field.  It  can  be  assumed  that  the  far 
field  is  still  governed  by  the  linearized  equation  for  unsteady 
flows,  only  the  assumption  that  the  waves  are  sinusoidal  in  time 
will  be  abandoned. 

Asymptotic  far  field  conditions  for  such  flows  are  obtained 
in  the  following  manner.  The  expression  Eq.(58')  combined  with 
Eq.  (6)  gives 

}  =  (f^sjp-1/2  +  f2(s)p"3/2)exp(iw(t-p)  (75) 

where  s  denotes  a  parameter  which  is  constant  along  a  characteristic 
and  p  identifies  the  station  along  the  individual  characteristics; 
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p  has  been  identified  with  the  local  radius  of  curvature  R  of 
the  lines  g  =  const.  At  the  moment,  it  is  assumed  that  one 
knows  the  angle  a  which  gives  the  direction  of  the  normal  to  the 
wave  fronts  and  therefore  also  the  angle  $  which  gives  the  direction 
of  the  characteristics  (Eq.  (64) ). (Notice  that  Eq.  (64)  does  not 
contain  the  frequency  uj)  .  We  denoted  by  d&  the  line  element 
along  a  characteristic.  We  have  found  in  Eqs.  (54)  and  (65) 

dp/dJl  =  dR/dJl  =  (M2  +  2M  cos  a  +  1)  ”1</2 

_i  (76) 

=  (M  cos  6  +  cos  y) 

For  simplicity,  we  omit  temporarily  in  Eq.  (75)  the  con¬ 
tributions  of  f 2*  Differentiating  Eq.  (75)  with  respect  to  t  and 
d&  one  obtains 

3<{>/3t  =  f^iojp  1//2exp(ico  (t-p) )  (77) 

3<J>/3&  =  fj.Op/34)  [-ia)p“1/2  -  (l/2)p"3/2]exp(io)(t-p) ) 

Combining  the  equations  with  Eq.  (75)  and  writing  R  instead  of  p, 
one  obtains  the  following  equation  which  is  free  of  u 

( 3«j»/32,)  +  (3R/3JI)  [  (3<J>/3t)  +  (2R)_1<{>]  =  0  (78) 

Equation  (78)  is  obtained  directly  from  Eq.  (61) ,  if  one  remembers 
that  each  factor  i«  arises  by  a  differentiation  with  respect  to  t. 
This  observation  allows  one  to  derive  far  field  conditions  which 
take  higher  order  terms  of  the  asymptotic  development  into  account 
immediately  from  Eq.  (60) 

[TT  +  IT  (2R  +  3t>  ]  [3l  +  “3l(  2R  +  aTt>  =0 

or  in  more  detail  from  Eq.  (63) 
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(79) 


32$/342  +  (3R/32.)  [(3/R)  (3$/34)  +  2  (32$/3*3t)] 

+  ( 3R/3A) 2 [ ( 3/4R2) $  +  (3/R) (3$/3t)  +  32$/3t2]  =  Q 

2  2 

Again  Eq.  (67)  can  be  used  to  express  3  $/3£  ,  which  in  the 

2  2 

notation  used  there  is  identical  with  3  $/3£  ,  in  terms  of 
derivatives  along  the  boundary  of  the  computed  part  of  the  flow 
field. 

So  far,  the  angle  a,  which  gives  the  direction  of  the  normal 
to  the  lines  g  =  const  has  been  assumed  to  be  known.  Here 
assumptions  analogous  to  those  of  the  preceding  section  must  be 
made,  if  one  assumes  that  the  origin  of  the  waves  is  known  (some 
point  in  the  vicinity  of  the  profile,  then  one  obtains  the  necessary 
formulae  by  replacing  in  Eqs.  (70),  (71),  and  (72),  ioi  by  3/3t. 

One  obtains 

( 3$/3Jl)  +  (g/R1)  (3$/3t)  =  0 

(3$/3i)  +  (g/R1) (3$/3t)  +  (1/2R1)$  =  0 

0Va*2)+  O/R^OS/aj l)  +  2(g/R1)  (32$/3£9t)  (80) 

+  (3/4rJ)$  +  (3g/R2)  3$/3t  +  (g/Rx)  2  (3  V  at2)  =  0 

It  is  possible  to  derive  in  this  general  case, at  least 
approximately,  the  normal  to  the  wave  fronts.  In  the  case  that  one 
has  sinusoidal  perturbations ,  we  assumed  that  along  a  line  g  = 
const  the  amplitude  of  <J>  and  also  its  phase  changes  only  slowly. 
Under  the  present  circumstances  this  amounts  to  the  assumption 
that 


grad  <|>  grad  g. 


One  then  obtains 


a  =  arctg  (<l>y/<J>x) 


(81) 
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This  defines  a  except  if  simultaneously  <J>  =  $  »  0.  Excluding 

y  x 

the  vicinity  of  such  points,  one  obtains  the  value  of  a  and 
therefore  also  the  local  values  of  p  in  the  same  manner  as  before. 

One  can  then  apply  directly  Equations  (74)  viz 

R-1  =  ay(l  +  M  cos  a)  (M  +  cos  a)-1  (82! 


R  1  =  -  <*x(.l  +  M  cos  a)  (sin  a)  ^ . 

Alternatively,  one  can  use  one  of  the  Eqs.  (73)  to  determine 
R.  The  factor  of  proportionality  between  g  and  <t>  cancels  under 
the  assumption  made  here;  one  can  replace,  for  instance, 


gxx/  I  grac*  g  I  bY  ^  I  g rad  <j>  |  . 


SECTION  VI 


CHARACTERISTIC  SURFACES  IN  THE  x,y,t  SPACE 

The  derivation  of  far  field  conditions  for  unsteady  flows 
shown  in  the  preceding  section  is  somewhat  indirect.  The 
following  discussion  is  an  attempt  to  throw  light  from  a  different 
direction  on  these  results .  As  far  as  the  practical  aspects  are 
concerned,  it  adds  nothing,  but  it  gives  the  connection  of  the 
present  approach  with  other  mathematical  concepts. 

One  observes  that  the  surfaces  of  the  x,y,t  space  given 
by 

t  =  g(x,y)  +  const  (83) 

are  characteristic  surfaces  for  the  partial  differential  equation 

(Eq.(5)).  The  derivation  is  found  in  Appendix  IV.  If  along  a 

noncharacteristic  surface  the  values  of  A,  $  ,  and  <j>.  are  given, 

x  y  u 

then  one  can  determine  by  means  of  the  differential  equation  the 
second  (and  higher)  derivatives  of  <J> ,  in  terms  of  derivatives 
formed  within  the  surface.  In  contrast,  these  derivatives  are  not 
uniquely  determined  for  a  characteristic  surface.  As  a  consequence 
it  is  possible  that  along  a  characteristic  surface,  discontinuities 
of  the  second  and  higher  derivatives  of  $  will  occur.  At  a  fixed 
point  (x,y,t)  the  orientations  of  these  surfaces  are  determined  by 
the  differential  equations  for  g,  namely 

(1-M2) g2  +  2Mgx  +  gy2  =  1  (84) 

We  introduce  coordinates  £,n  within  one  of  the  characteristic 
surfaces,  £  and  g  are  identical  with  the  x  and  y  coordinates  of 
a  point  (x,y,t)  of  the  surface.  Let 

u  ■  *x,  V  =  ♦  ,  W  =  <frt.  (85) 
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Then  one  has,  (since  t  =  g(x,y)) 


^xx  +  ^xt^x 

vn  ■  ♦yy  +  *ytgy  <8S 

w5  “  *xt  +  *tt9x 

W  =  *  .  +  d> ,  .  g  . 
n  vyt  vtt’y 

In  order  for  second  derivatives  of  $  to  exist,  the  values  of  U,  V, 
and  W  within  the  surface  must  be  connected  by  the  compatibility 
condition 


-(1-M2)U^  -  +  (2M  +  (1-M2)  gx)  W  ^  +  gyW^  =  0.  (87) 

This  derivation  is  found  in  Appendix  IV. 

A  plot  of  lines  g  =  const  in  the  x,y-plane  can  be  interpreted 
in  two  different  manners.  If  one  considers  t  in  Eq.  (83)  as  fixed, 
then  these  curves  give  the  intersection  of  characteristic  surfaces 
in  the  x,y,t  space  with  a  plane  of  constant  t.  If  one  considers  the 
constant  in  Eq.  (83)  as  fixed,  then  the  curves  represent  a  contour 
map  of  one  characteristic  surface  where  t  corresponds  to  the 
altitude . 

Assume  now  that  for  a  given  time  a  curve  g(x,y)  =  cQ  is 
given.  To  construct  the  characteristic  surface  one  must  find  the 
curves  g  =  const  for  other  values  of  t.  For  this  purpose  one  can 
use  the  results  found  above.  Consider  points  (x,y)  reached  from  a 
certain  point  (xQ,yo)  of  the  starting  curve  g(x,y).  =  cQ  by  traveling 
along  straight  lines  given  by 

(x-xQ)  =  (p-pQ)(M  +  cos  a) 

(y-yQ)  =  (P-P0)sin  a 
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where  a  is  the  angle  of  the  normal  to  the  curve  with  the  x-axis. 

These  are  the  directions  given  by  Eq.  (52) . 

Since,  according  to  Eq.  (54) 

dg/dp  =1  (89) 

one  obtains  a  point  of  the  curve  g  =  c^  by  setting  p-pQ  =  c^-cQ. 

The  same  point  is  reached,  if  one  proceeds  from  the  original 

point  in  the  x-direction  by  a  distance  M(c^-cq)  and  then  in  the 

direction  of  the  normal  by  the  distance  c.-c  .  This  is  in 

l  o 

accordance  with  the  description,  given  in  Section  II,  that  a  wave 
front  spreads  as  in  still  air  except  that  the  perturbation  floats 
downstream  with  the  (dimensionless)  velocity  M. 

For  a  further  discussion,  we  consider  all  possible  values  of 
gx  and  gy  (compatible  with  Eq.  (84)),  they  are  obtained  by  letting 
a  vary  in  Eq.  (88)  from  0  to  2tt.  The  locus  of  the  points  g  =  ^ 
which  are  reached  from  point  xQ,yo  is  then  given  by 

x  =  x^  +  (c.  -  c  )  (M  +  cosa) 

°  1  °  (90) 

y  -  yQ  +  (cx  -  cQ) sin  a. 

This  is  a  circle  with  center  x  +  M(c,  -  c  )  and  radius  (c.  -  c  ) . 

o  1  o  1  o 

In  the  x,y,t  space,  these  circles  are  the  contour  lines  of  an 
oblique  cone  with  tip  at  point  ( xQ,yo )  and  a  generator  given  by  the 
above  circle.  At  the  point  of  the  pertinent  characteristic 

x  =  x^  +  (c-,  -  c  )  (M  +  cosa) 
o  x  o 

and  y  =  y  +  (c.  -  c  )sin  a 

o  l  o 

the  circle  drawn  for  a  fixed  point  of  the  line  g  =  c^  and  the  line 
g  =  c^  have  the  same  tangent;  consequently,  the  cone  mentioned 
above  is  tangent  to  the  characteristic  surface  in  the  x,y,t  space. 

The  lines  of  tangency  are  called  bi-characteristics.  The  characteristics 
in  the  x,y-plane  with  which  our  considerations  started  are  the 
projections  into  the  x,y-plane  of  the  bi-characteristics  in  the 
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x,y,t  plane.  The  line  g  =  c  can  he  obtained  by  letting  the  initial 
line  float  in  the  x-direction  by  the  distance  M(c^  -  cQ) ,  by  drawing 
circles  with  radius  c^  -  cQ  and  centers  at  the  curve  so  shifted,  and 
finally,  by  determining  the  envelope  of  these  circles.  According  to 
the  description  given  here,  each  characteristic  surface  is  determined 
by  the  initial  line  g  =  const.  If  two  such  lines  are  tangent  to 
each  other  at  some  point  (xQ,yo) /  then  they  will  be  tangent  to  each 
other  along  the  bi-characteristic  that  starts  at  this  point. 

It  is  natural  to  introduce  in  the  compatibility  condition 

(Eq.  (87))  instead  of  £  and  n,  the  derivatives  in  the  direction 

of  the  bi-characteristics  and  of  the  curve  g  =  c  .  We  denote  by 

o 

£• the  projection  of  the  bi-characteristic  into  the  x,y-plane, 
and  by  r^,  the  coordinate  in  the  direction  of  the  tangent  to  the 
line  g  =  const.  The  projection  of  the  bi-characteristic  into  the 
X/Y-plane  is  identical  with  the  characteristic  lines  considered 
above.  The  angle  of  with  the  £-axis  (which  coincides  with  the 
x-axis)  is,  therefore,  given  by  the  angle  8,  determined  by 
Eq.  (43) .  The  angle  of  the  normal  to  the  contour  is  denoted  as 
before  by  a.  One  then  has 


£  = 

£^cos  6  -  sin  a 

n  = 

£1sin  6  +  n1  cos  a. 

3 

3*1 

=  cos 

r 

+  sin  8~!? 

3 

=  sin 

3 

aTT 

J  3 

+  cos  a  — t— 

3rj 

and 


(91) 


(92) 
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3  3 

(cosaj| - sinB-^pj— ) /cos  (8-a) 

Csinaj|-  +  cosBypp) /cos  (B-a)  . 


(93) 


This  and  the  expression  for  gx  and  g  from  Eq.  (50)  are  now 
substituted  into  Eq.  (87) ;  one  obtains 


-(1-fT)  cos  aOU/3^)  -  sina(3V/3C1)  +  (1  +  M  cosa)  OW/3^) 

+  (M2+2M  cosa  +  1)  ”1,/2[sina(l-M2)  (3U/3n1)-(M  +  cosa)  OV/Sr^) 


-M  sina ( 3W/3n1) ]  =  0 


(94) 


If  one  changes  a  by  ir  then  one  obtains  a  second  bi¬ 
characteristic.  Let  £2  be  its  direction.  The  corresponding  value 
of  8  (denoted  by  82)  is  then  given  by 


2  -1/2 

cos  32®  (M  -  cos  a) (1  -  2M  cos  a  +  M  ) 

2  -1/2 

sin  82  =  “  sina(l  -  2Mcos  a+  M  )  '  . 


(95) 


The  direction  dn1  is  changed  into  -dri^.  One  then  obtains 


(1-M  )cosa(3U/3?2)  +  sina(3V/3C2)  +  (1-Mcos  a)(3W/3?2) 

+  (M2-2Mcos  a+  l)~1/2[sina(l-M2)  (SU/Sn^  (M-cos  a)(3V/3n1) 


(96) 


-  M  sinaOW/Sn^l  -  0. 


In  preparation  for  a  discussion  of  far  field  conditions,  we 
now  discuss  the  following  configuration.  We  consider  a  characteristic 
surface  t  =  g(x,y)+cQ  (with  cQ  fixed),  in  the  vicinity  of  a  point 
x,y.  For  time  t  >  g(x,y)+  cQ  let  $  be  identical  to  zero.  At  the 
characteristic  surface  as  everywhere  else,  $  is  continuous,  but  its 
first  derivative  changes  very  rapidly  as  far  as  this  is  compatible 
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with  continuous  <J>.  For  t  >  g(x,y)  +  cq  +  £  ,  the  first  derivatives 
reach  some  fixed  (or  nearly  fixed)  value.  Ultimately,  we  allow  £ 
to  tend  to  zero,  then  a  jump  of  some  first  derivative  propagates 
along  the  characteristic  surface.  In  this  situation  the  derivatives 
of  U,  V,  and  W  (4>x,  and  4>t)  with  respect  to  n  (that  is  within 

the  characteristic  surface) ,  remain  bounded.  The  derivatives  along 
the  bicharacteristic  which  belongs  to  the  same  normal  to  the 
characteristic  surface  and  which  passes  through  this  surface 
becomes  very  large.  We  integrate  Eq.  (96)  with  respect  to  £2, 
the  path  of  integration  goes  from  the  undisturbed  region  through 
the  jump  region  to  a  point  shortly  beyond  it.  Let  [U] ,  [V] ,  and 
[W]  denote  the  jumps  of  these  quantities;  they  are  actually  identical 
with  U,  V,  and  W  at  the  point  beyond  the  jump  region.  Let  the 
expression  containing  derivatives  with  respect  to  ri  temporarily 
be  denoted  by  Q2«  Q2  remains  bounded  even  in  the  jump  region. 

One  then  obtains 

(1-M2)  cosa  [U]  +  sinci[V]  +  (1-M  cosa)  [W]  +  /Q2d£2  =  0.  (97) 

In  the  limit  as  i  -*  0 ,  the  length  of  £2  over  which  the  integration 
is  extended  becomes  zero.  One  thus  obtains 

(1-M2)  cosa*  +  sinacj)  +  (1-M  cosa)  <f>.  =  0.  (98) 

x  y  z. 

Because  of  the  continuity  of  <j>,  and  since  <j>  ■  0,  before  the  jump 
the  gradient  of  <p  is  normal  to  the  curve  g(x,y)  =  t-cQ,  for  any 
fixed  time.  One  remembers  that  a  gives  the  normal  to  lines  g(x,y)= 
const .  Therefore , 


<b  =  |grad<J>|  cos  a 

x  (99) 

4>y  ■  |g3rad<|>|  sin  a. 

Thus,  one  obtains  from  Eq.  (98) 

(1  +  Mcosa)  |grad<f>J  +  4>t  *  0.  (100) 
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We  introduce  the  direction  of  the  characteristics  in  the  x,y- 
plane;  we  have  found  that  such  characteristics  are  the  projection 
of  the  bicharacteristics  (of  the  characteristic  surface)  into  the 
x,y  plane.  The  angle  of  the  characteristics  with  the  x-axis  has 
been  denoted  by  0;  the  line  element  of  the  characteristic  by  d£. 

Since  <j>  =  const  along  the  line  g(x,y)  =  t-cQ,  one  has 

=  |grad<j)|  cos(0-a)  (101) 

and  with  Eq.  (64) 

||  =  |  grad<j>  j  ( 1+M  cos  a)  (M2+2M  cos  a+  D  "1/2  (102) 

Equation  (100)  can  then  be  written  in  the  form 

(3<j>/3£)  +  .  (M2  +  2M  cosa  +  1)  "1/2(3<J>/9t)  =  0  (103) 

This  is  Eq.  (78)  combined  with  Eq.  (77)  but  with  the  term  (2R  1) 
omitted.  For  a  jump  of  (grad<f>|  along  a  characteristic  surface 
the  last  formula  is  exact. 

In  general,  the  expression  / Rd^  which  vanishes  for  a  jump 
will  play  a  role.  The  line  along  which  one  integrates  intersects 
the  characteristic  surfaces  t  =  g(x,y)+c.  In  an  oscillatory  flow 
these  are  surfaces  of  constant  phase.  But  for  u>  large  the  phase 
changes  rapidly  as  one  passes  from  one  such  surface  to  the  next 
one  (this  happens  if  one  moves  in  the  direction  ) •  T^e  derivatives 
with  respect  to  (which  are  derivatives  along  the  different  line 
g  *  const)  remain  bounded.  One  thus  obtains  an  integral  with  a 
strongly  oscillatory  integrand;  as  w  increases  the  oscillations 
become  narrower.  Hence  it  follows,  that  in  a  low  order  approxima¬ 
tion  in  uf1,  the  integral  can  be  disregarded.  It  is  perhaps 
possible  to  estimate  this  integral  (using  integration  by  parts 
and  the  second  mean  value  theorem) ,  but  actually  the  derivation 
of  higher  order  formulae  in  the  manner  shown  in  Section  V  is 
simpler. 


The  last  discussions  have  shown  that  the  far  field  conditions 
can  be  interpreted  as  an  expression  for  the  compatibility  condition 
for  the  bi-characteristics  of  a  special  characteristic  surface, 
defined  by  the  postulate,  that  for  a  fixed  t  it  intersects  the 
characteristic  surface  of  constant  phase  t  *  g(x,y)  +  c  ,  along 
the  curve  g(x,y)  =  t  -  cQ.  This  special  choice  is  necessary  in 
order  to  ensure  that  the  factor  R  in  the  integral  / Rd^2  remains 
bounded. 

In  Eq.  (103)  the  radius  of  curvature  of  the  line  g(x,y)  =  const 
does  not  appear.  The  results  agree  with  Eq.  (78)  only  for  R  •+**». 

In  other  words,  one  uses  a  relation  for  plane  waves.  The  dependence 
on  the  parameter  p  (which  is  identical  with  the  radius  of  curvature 
R)  and  which  would  be  needed  in  order  to  derive  higher  order 
approximations  can  be  obtained  on  the  basis  of  the  last  formula, 
if  one  uses  in  addition  the  compatibility  conditions  for  the 
surface  t  =  g(x,y)  +  cQ. 

First  we  express  the  fact  that  at  this  surface  4>  =  const. 

Since  a  gives  the  direction  of  normal  to  the  line  g(x,y)  in  the 
x,y-plane  and  since  4>  =  const,  one  has: 


<j>x  =  |  grad<J)  |  cosa 
<|>y  =  |  grad<J)  |  sina 


(104) 


where  |grad<ft|  refers  to  the  function  <j>(x,y,t)  for  constant  t. 
Along  the  intersection  of  a  plane  y  =  const  with  t  =  g(x,y)  +  cq 
one  has 


dt  =  gxdx. 


Therefore,  since  <p  m  0 


d>,  +  <b  3?  =  <p  +  g  ^6  . 

Tt  yx  at  Tt  yx  Yx 


(105) 
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F* 


Hence , 


|grad<J>|cos  a 


and  with  Eq.  (50) 


4>t  =  -  j  gradtf)  |  (1+m  cosa) 


(106) 


Our  next  task  is  the  determination  of  the  derivatives  with  respect 
to  r\^  in  Eq.  (94)  .  The  expression  to  be  evaluated  is 

Ql  =  [sina(l~M2) (3U/3n1)  -  (M+cosa) ( SV/Sr^)  -  Msina ( SW/Sr^) ) 

[M2  -  2Mcosa  +  1]1/2  (107) 

with 


U  =  =  j  grad<J>  |  cos  a 

v  =  <J>y  =  |  grad<j»  j  sin  a 
W  =  <pt  =  |  grad<j>  |  ( 1+M  cosa). 


The  derivatives  of  U,  V,  and  W  in  the  expression  Q^,  therefore,  are 
expressed  in  terms  of  derivatives  of  j  grad<J>  j  and  of  a  with  respect 
to  r\1  The  coefficient  of 


|  grad*  | 

turns  out  to  be  zero,  for 
2 

sina(l-M  )cosa  -  (M  +  cosa) sina  +  Msina (1+Mcosa)  =  0, 


Therefore , 


Ql  •  -jgrad^l^— “  ( 1+Mcosa)  (M2  +  2Mcosa  +  1)1/2. 
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In  a  similar  manner,  one  obtains  for  the  derivatives  in  Eq.  (94) 
-(l-M2)cosaOU/H1)  -  sinaOV/9^)  +  (1+Mcosct)  ( SW/3^)  = 

-2  [3  |  grad^>  | /3C1J  (1+Mcosa) 

No  terms  3a/3£^  occur.  This  is  in  keeping  with  Eq.  (32) ,  which 
shows  that  along  a  characteristic  (here  bi-characteristic)  gx,  g^ 
and  consequently  also  a  are  constant.  One  also  observes  that 

da/dri  ^  =  R  1 

(where  as  before  R  is  the  radius  of  curvature  of  a  line  g  =*  const)  . 
One  thus  obtains  from  the  compatibility  condition  Eq.  (94) 

-2(d|grad4>|/d£1)-R_1|grad(})  j  (M2  +  2Mcosa  +  1)"1/2  «  0 

In  Appendix  II  we  have  identified  R  with  a  parameter  p, 
the  line  element  d£^  corresponds  to  d£.  Using  Eq.  (53)  one  then 
obtains 

d  ( log  |  gradij)  | ) /dp  >  (l/2p)  *  0. 


Hence , 

|grad<J>|*  const  p-*^2  *  const  R  ^'/2 

Prom  Eqs.  (104)  and  (105) 

*  -1/2 

<j>x  =  const  R  '  cos  a 

-1/2  . 

*  const  R  sm  a 
■  const  R  ^2(1  +  Mcosa) 


(108) 


(109) 
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The  constant  in  these  three  expressions  is  the  same  but  may 
vary  from  characteristic  to  characteristic.  According  to  Eq.  (77) 

dR/di.  =  (M2  +  2Mcosa  +  1)“1/2. 

Note  that  these  expressions  satisfy  the  Sommerfeldt  radiation 
condition  (Eq.  (78) )  with  the  term  C2R  ^)  omitted.  This  can  be 
shown  with  the  aid  of  Eqs.  (102)  and  (77) . 

The  conditions  (103)  arise  from  bi-characteristics  which 
cross  the  characteristic  surface  along  which  a  jump  of  |grad<j>| 
propagates.  Here  the  contribution  of  Q2  vanishes  in  the  limit 
of  a  sudden  change  of  |grad<j>|.  In  this  equation,  no  contributions 
of  R  appears.  The  expressions  Eq.  (108)  and  (109)  are  derived 
from  the  relations  for  bi-characteristics  lying  directly  above 
the  characteristic  surface  at  which  the  jump  occurs.  Here  the 
term  Q^,  the  counterpart  of  Q^,  is  not  negligible  in  comparison 
to  other  terms  of  the  compatibility  condition.  This  fact  is 
responsible  for  the  dependence  upon  the  radius  of  curvature  R 
of  the  line  g(x,y)  =  const. 
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SECTION  VII 

THREE-DIMENSIONAL  PROBLEM 


The  three-dimensional  problem  has  actually  more  practical 
significance  than  the  two-dimensional  problem.  Most  of  the 
necessary  formulae  can  be  fairly  easily  derived  by  using  the 
analogy  with  the  two-dimensional  problem.  Eq.  (22)  is  now  replaced 
by 


<p(x,y,z)  =  h(x,yz)  exp  (-iwg  (x,y ,  z) )  .  (110) 

Equation  (24)  is  replaced  by 

(1-M2)g2  +  g2  +  g2  +  2Mgx  -1=0.  (Ill) 

Let  D  denote  differentiation  in  the  direction  of  a  characteristic,  and 
let  D£  be  the  line  element  of  a  characteristic.  With  unit  vectors 
ex,  ey,  and  ez  respectively  in  the  x,y,  and  z  directions,  one 
obtains 


e  (Dx/Dl)  +  e  ( Dy/DZ )  +  e  (Dz/DJL)  = 
x  y  z 

=  const  [ex(  (1-M2)  gx+M)  +ey9y  +  ezg2J 
Then  one  finds,  that  along  the  characteristics 


(112) 


Dgx/Di  =  Dgy/Di  =  Dgz/DA  =  0. 

Therefore,  gx#gy,  and  gz  are  constants.  The  characteristics  are 
straight  lines.  Let  (n,x) ,  (n,y) ,  and  (n,z)  be,  respectively,  the 
angles  of  the  normal  to  the  surface  g  =  const  with  the  x,  y,  and  z 
axes.  One  has,  of  course, 

cos2(n,x)  +  cos2(n,y)  +  cos2(n,z)  =  1  (113) 
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Equation  (50)  suggests  that 


jgrad  g|  =  (1  +  Mcos(n,x))~ 
gx  =  cos(n,x)(l  +  Mcos(n,x)) 
gy  =  cos(n,y)(l  +  Mcos(n,x)) 
g_  =  cos(n,z)(l  +  Mcos(n,x)) 

Z 


(114) 


On  readily  verifies  that  this  is  the  solution  of  Eq.  (Ill) . 

It  depends  upon  the  parameters  cos(n,x),  cos(n,y),  and  cos(n,z). 
Because  of  Eq.  (113),  there  are  actually  only  two  free  parameters. 
One  also  verifies  that  Eq.  (112)  is  satisfied  if  one  sets 

x  *  xq  +  p  (cos(n,x)  +  M) 

y  =  yQ  +  P  cos (n,y)  (115) 

z  =  z  +  p  cos(n,z) 
o 

where  p  is  a  parameter  whose  meaning  so  far  is  not  specified  (in 
the  two-dimensional  case,  p  has  been  identified  with  the  radius  of 
curvature  of  the  curve  g  =  const  at  the  point  under  consideration; 
in  the  three-dimensional  case,  one  does  not  have  a  single  radius  of 
curvature) .  One  then  has 

DJI/Dp  =  (M2  +  2Mcos  (n,x)  +  1)1//2  (116) 

Dg/Dp  =  gx(cos(n,x)  +  M)  +  g^cos(n,y)  +  gzcos(n,z)  =  1  (117) 

and  hence 

(Dg/Di)  =  (Dg/Dp)  (Dp/D 2.)  *  (M2  +  2Mcos(n,x)  +  1)"1/2. 

If  all  the  characteristics  emanate  from  the  origin,  then 
one  has 


(118) 


x  =  p  (cos(n,x)  +  M) 
~y  =  p  cos (n,y) 
z  =  p  cos (n,  z) ; 


(119) 


g  =  gvx  +  g  y  +  g_z  =  p  .  (120) 

a  y  z 

A  surface  g  =  const  is  a  sphere  with  radius  p  and  center  M  p,  for 
one  has  as  a  consequence  of  Eqs.  (114)  and  (119) 

(x  -  Mp)2  +  y2  +  z2)  =  p2.  (121) 

In  analogy  with  Eq.  (36) ,  we  set , tentatively / for  this  case 

g(x,y,z)  = (1-M2)"1[-Mx  +  (x2  +  (1-M2) (y2  +  z2))1/2].  (122) 

Then 

gx  =  (1-M2)_1{-M  +  x(x2  +  (1-M2) (y2  +  z2) ) ~1/2] } 
gy  *  y [x2  +  (1-M2) (y2  +  Z2)3'1/2 
gz  -  z[x2  +  (1-M2)  (y2  +  z2)]'1/2. 

Then  one  can  verify  that  Eq.  (Ill)  is  satisfied.  Next  one 
must  determine  H.  In  principle,  one  could  return  to  the  counter¬ 
part  of  Eq.  (26)  and  evaluate  the  coefficients  of  this  equation. 

But  one  can  also  refer  to  the  counterpart  of  the  particular 
solutions,  Eq.  (16).  They  are  spherical  harmonics,  and  they 
have  the  form 

r  1P(r”1)exp(-iyr) 

with 

r  =  (x2  +  y2  +  z2)1/2.  (123) 
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P  is  the  generic  expression  for  a  power  series.  The  analog  to 
Eq.  (58)  suggests  the  following  approximations  for  <j> 

(124) 

(J)  =  If1(x/R1,y/R1/z/R1)  r"1  +  f2(x/R1,y/R1,z/R1)  R^2  ]  exp  (-iog  (x,y ,  z)  ). 

Along  a  characteristic , the  functions  f^  and  f2;along  with  the 

values  of  g  ,g  ,  and  g  .are  constant.  This  leads  to  the  counter- 

■X  y  z 

part  of  Eq.  (60) 

*57  +  R^  +  iw  51^  *57  +  R^  +  iaj  D?}  0  (125) 

and  the  simpler  (and  less  accurate  forms) 

[Dl  +  +  -§fH  “  0  (126) 

and 

Iq?  +  io>  gf-]*  =  o.  (127) 

Here 

Rx  =  (x2  +  y2  +  z2)1/2  (128) 

The  direction  of  the  characteristics  is  given  by  the  rays  through 
the  origin.  One  has 


22.  .  _a 

DZ  Rx 


and  g  is  given  in  Eq.  (122) . 


(129) 


For  general  unsteady  flows  with  waves  originating  from  the 
origin,  one  than  obtains  the  far  field  conditions 


_3 

3t 


+ 


+ 


0.- 


(130) 
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or  the  simpler  forms 


Hjt  +  -^+  =°  <131> 

'-§r  +  ‘0-  <132> 

The  second  derivative  which  occurs  in  the  above  equation  is 
expressed  by 

D2/DZ2  =  (x2/R2) (32/3x2)  +  (y2/R2) (32/3y2)  +  (z2/R2) (32/32) 

111  (133) 

+2(xy/R^) (32/3x3y)  +  2 (xz/R2) ( 32/3x3z)  +  2 (yz/R2) ( 32/3y 3z) . 

Assume  that  one  boundary  of  the  computed  flow  field  is  a  plane 

x  =  const.  From  the  values  of  $  and  <J>x  within  this  plane,  one 

determines  <j>yy,  <j>22,  *  ,  <t>xz,  4>xy,  and  ^  is  expressed 

in  terms  of  <j>  ,  <f>  ,  $  .  ,  and  <p  by  the  differential  equation 

yy  zz  At  1 1 

for  <(). 

Equations  (124)  through  (133)  give  far  field  conditions  for 

which  one  estimates  the  origin  of  the  waves  (in  essence,  the 

surfaces  g  =  const) .  The  location  of  the  origin  determines  the 

direction  of  the  characteristics  (the  ray  through  the  origin  of  the 

222  1/2 

wave) ,  the  value  of  R^  =  (x  +  y  +  z  )  '  .  One  obtains 

g/Rx  =  (l-M2)"1!-*!*  +  (x2+(l-M2)  (y2+z2)1/2]  [x2+y2+z2] -1/2 . 

This  approximation  may  sometimes  be  unsatisfactory.  In  three- 
dimensional  problems  at  an  intermediate  distance  from  the  airplane, 
surfaces  g  =  const  are  more  likely  to  resemhle  ellipsoids  than 
spheres.  Then  one  may  proceed  as  follows.  We  consider  sinusoidal 
perturbations.  One  identifies  surfaces  g  =  const  with  surfaces  of 
equal  phase  of  <J>.  One  determines  by  numerical  differentiation 
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(134) 


cos(n,x)  ■  | grad  g| 

cos  (n,y)  =  g  / )  grad  g}- 
cos (n, z)  =  g  /| grad  gj . 

One  ought  to  have,  according  to  Eq.  (114) 

|grad  gj  =  (1  +  Mcos(n,x))-^  (135) 

This  gives  a  check  (and  probably  also  an  adjustment)  of  the  numerical 
results.  The  characteristic  direction  then  follows  from  Eq.  (115) 

Dx/Di  **  (cos(nx)  +  M)/(M2  +  2Mcos(n,x)  +  1)  ^2 
Dy/DS,  *  cos(ny)/  (M2  +  2Mcos(n,x)  +  1)  ^2 
Dz/DJl  *  cos(nz)/  (M2  +  2Mcos(n,x)  +  1) 1//2 

This  defines  all  quantities  occurring  in  the  Sommerfeldt 
far  field  condition,  (the  generalization  of  Eq.  (132)) 

D<j>/D&  +  iwDg/DJl$  «  0. 


The  next  higher  approximation  which  corresponds  to  Eq.  (131)  requires 
the  evaluation  of  the  function  h  in  Equation  (110) .  Proceeding  in 
analogy  with  the  two-dimensional  case  (Eqs.  (25)  and  (26)),  one 
arrives  at  the  analog  to  Eq.  (34) 


1/2 

DdogH^/DA  +  (1/2)  [l-M2)gxx+gyy+gzz]/[g2+g2  +((1-M2)gx+M)2]  -0 . 

Using  Eqs.  (134)  and  (135) ,  one  obtains 


(136) 

D(logH.)/D£  +(1/2) [1-M2) g  +g  +g  ] [l+Mcos(nx) [M2+2Mcos (nx) +1] "1/2 . 

X  XX  yjf  Z  Z 


Then  one  obtains  from  the  approximate  representation  for  $ 
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<p  =  H1  exp(-itjg) 

D<p/DZ  +  [iw  (Dg/DZ)  -  D  (logH1)  /DZ]  <j>  =  0. 

(logH.)D  is  evaluated  from  Eq.  (.136) .  This  requires  the  determination 

2 

of  the  factor  (1-M  }g  +  g  +  g  which  can  be  done  numerically. 

A  geometrical  interpretation  is  formed  in  Eq.  (A. 31),  but  it  is 

preferable  to  express  directly  that  along  the  characteristics,  gx, 

g  ,  and  g  are  constant. 

’y  ^z 

gxx(M  +  cos(n'x))  +  gxycos(n,y)  +  gxzcos(n,z)  =  0 

gxy(M  +  cos (n,x) )  +  gyycos(n,y)  +  gyzcos(n,z)  =  0 

g„_  (M  +  cos  (n,x) )  +  g  cos(n,y)  +  g„cos(n,z)  =  0. 

az  yz  zz 

We  describe  the  procedure  for  a  boundary  surface  given  by  a  plane 

x  =  const.  There  one  finds  g  ,  g  ,  and  g  by  numerical 

x  y  Z 

differentiation.  This  allows  one  to  determine  cos(n,x),  cos(n,y), 
and  cos(n,z) .  By  further  differentiation  within  this  plane  (that 
is  with  respect  to  y  and  z)  one  can  express 


Jyy' 

Then  one  obtains  g%„.  and  g 

xy  ’xz 

gxy  -  -(gyy  cos (n,y)  + 
gX2  *  -<9yz  cos (n,y)  + 
and  from  the  first  equation 


g  ,  and  g 
’yz'  ’zz 

from  the  last  two  equations 
gyzcos(n,z) )/(M  +  cos(n,x)) 
g„,cos(n,z) )/(M  +  cos(n,x)) 

Z  Z 


2 

gxx  “  ^yycos  (n'y)  +  2gyzcos(n,y)cos(n,z) 
This  allows  the  evaluation  of  d(logH^)/d£, 


(137) 

2 

+  g__cos  (n,z)/M+cos(n,x) ) 
z  z 

Eq.  (136) . 


2 
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An  interpretation  of  the  expression 


(1/2)  tl”M2)gxx  +  gyy  +  g2zJ [1  +  Mcos(nx)] 

which  occurs  on  the  right  hand  side  of  the  equation  for  DdogH^DH 
is  given  at  the  end  of  Appendix  II.  Using  Eq.  (135)  one  recognizes 
that  it  represents  the  average  curvature  of  the  surface  g  =  const 
at  the  point  under  consideration.  The  average  curvature  is 
obtained  in  the  following  manner.  One  forms  the  curvature  of  two 
curves  which  arise  from  the  intersection  of  the  surface  g  =  const 
with  two  planes  through  the  normal  to  this  surface  which  are 
perpendicular  to  each  other,  and  then  forms  the  average.  One  of 
these  planes  can  be  chosen  arbitrarily,  the  other  one  is  then 
determined.  The  average  curvature  is  independent  of  the  orienta¬ 
tion  of  these  planes.  For  a  certain  orientation,  the  radii  of 
curvature  assumes  simultaneously  extreme  values .  Let  these  values 
be  Q  and  R2  ^  for  the  surface  g  =  gQ,  and  assume  that  for  this 
surface  p  =  0.  The  dependence  upon  the  parameter  p  along  a  fixed 
characteristic  is  then  given  by 

c average  ’  U/2>  CR10+P>  '^(Kjo+P)  _1- 

One  recognizes  that  the  expression  (136)  contains  certain  terms  of 
higher  order  in  p-1. 

With  Eq.  (116)  one  can  now  rewrite  Eq.  (136) 

D(logH1) /dp  *  -  |[(R10+p)“1+(R20+p)"1] . 


Hence , 


Hx  =  const  [ (R10  +  p) <R20  +  P)]"1/2.  (138) 

The  expression  HR^q+p)  O^q+P)  3  is  the  Gaussian  curvature  of 
the  surface  g  =  const  at  the  point  under  consideration.  The  constant 
may  change  from  characteristic  to  characteristic.  For  p  large  one 
obtains 


=  const  p 

This  is  the  form  suggested  by  the  particular  solutions  (123) . 

The  present  form  of  the  far  field  conditions  is  somewhat  more 
general.  In  the  two-dimensional  case  one  of  the  principal  radii 
of  curvature  (and  the  constant  in  Eq.  (138))  are  infinite.  Then 
one  obtains  the  results  of  the  preceding  section 

-1/2 

=  const  p  '  . 

The  estimation  of  the  H2  carried  out  in  the  two-dimensional 
case,  is  valid  only  for  p  large.  Then  one  has  for  the  three- 
dimensional  case  ^  R2,  and 

H2  v const  (riR2)-1- 

On  this  basis  it  would  be  possible  to  formulate  approximations  to 
the  far  field  conditions  of  the  next  order.  However,  one  may 
have  doubts  whether  in  a  practical  case  this  refinement  is 
justified.  The  derivations  are  based  on  the  assumption  that 
surfaces  of  constant  g  can  be  approximated  by  surfaces  of  constant 
phase  of  <p .  This  is  only  an  approximation.  The  application  of 
these  ideas  to  general  unsteady  flows  presupposes  again,  that  the 
computed  flow  permits  one  to  recognize  wave  fronts  (which  define 
the  surfaces  g  =  const) .  In  critical  cases  it  is  probably 
preferable  to  use  far  field  conditions  which  are  more  complicated 
but  do  not  require  the  assumption  of  high  frequency  or  the 
identification  of  wave  fronts  or  wave  origins. 
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SECTION  VIII 


OSCILLATORY  PERTURBATIONS  IN  A  PLANE  FLOW  FIELD  WITH 
A  FREE  STREAM  MACH  NUMBER  ONE 


Figure  2  shows  the  wave  patterns  obtained  in  the  linearized 
treatment  of  subsonic  flows,  which  arises  in  the  limit  where  the 
free  stream  Mach  number  goes  to  one.  The  perturbations  pile  up  at 
the  value  of  x  where  the  waves  originate.  This  happens  because  the 
sound  velocity  which  governs  the  manner  in  which  perturbations 
spread  through  the  flow  field  is  considered  as  constant  throughout 
the  flow  field.  In  reality,  the  profile  generates,  even  at  a  free 
stream  Mach  number  one,  a  subsonic  region  upstream  of  the  profile 
in  which  perturbations  will  travel  upstream.  A  pile-up  of  waves 
will  not  occur,  although  at  a  great  distance  from  the  profile  the 
wave  length  pertaining  to  a  certain  frequency  will  be  very  short. 

The  following  analysis  takes  the  modification  of  the  sound  velocity 
by  the  presence  of  the  profile  into  account.  At  a  sufficient 
distance  from  the  body  the  flow  field  at  a  free  stream  Mach  number 
one  can  be  described  by  a  similarity  solution.  In  the  following 
we  study  periodic  perturbation  in  this  part  of  the  flow  field. 

The  differential  equations  for  the  underlying  steady  flow 
field  is  obtained  in  the  familiar  manner;  one  assumes  that  the 
deviations  from  the  sonic  parallel  flow  are  small,  but  retains 
a  critical  term  which  allows  the  change  of  the  type  of  the 
differential  equation  from  elliptic  in  the  subsonic  region  to 
hyperbolic  in  supersonic  region.  The  simplified  differential 
equation  for  the  steady  flow  field  is  then  given  by 

--X&  4  <*ox,+  ♦oyy-  °-  (139> 

Here  Is  the  potential  which  gives  in  the  basic  steady  flow  field 
the  deviation  of  the  flow  from  a  parallel  flow  with  the  Mach  number 
one.  In  this  field  one  now  studies  unsteady  perturbations  which  are 
small  even  in  comparison  to  the  steady  perturbations  introduced  by 
the  presence  of  the  body.  Let  #  be  the  potential  describing  these 
perturbations.  They  satisfy  the  following  partial  differential 
equation 
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-<*+  +  V  -  2*xt  -  ht  “  °-  U40) 

Here  the  space  coordinates  have  been  made  dimensionless  by  a 
characteristic  length  L  and  the  time  by  (x,a)_1  (where  "a"  is  the 
velocity  of  sound.  The  basic  potential  is  considered  as  known. 
Eq.  (140)  is  the  counterpart  of  Eq.  (5) .  Introducing 

$(x,y,t)  =  <j>(x,y)exp(io>t)  (141) 


one  obtains 

•  -  (y+  l)-r-^(<j>.  <b  )  d>  +  d>  -  2iw$  +  u2(j)  =  0  (142) 

1  3x  yoxyx '  Yx  Yyy  yx  T 

which  corresponds  to  Eq.  (9) .  In  analogy  to  Eq.  (22)  one  sets 

4>  =  h(x,y)exp(-iuig(x,y) )  .  (143) 

Substituting  this  equation  into  Eq.  (142)  and  collecting  the  terms 
2 

with  powers  w  and  ui  one  obtains  the  counterparts  to  Eqs.  (24)  and 
(26) 


(y  +  u*oxgx  -  2gx  -  «y  +  1  =  0 


(144) 


and 


(145) 


(<Y+l)4ox2gx-2]<VH)  -  2gsr(Hy/H)  +  (Y+l)(4oxgxx++oxxgx)-gyy>0. 


As  in  Section  II  we  reduce  the  integration  of  Eq.  (145)  to  the 
integration  of  ordinary  differential  equations  by  introducing  the 
concept  of  characteristics.  Differentiating  Eq.  (144)  with  respect 
to  x  and  y  one  obtains 


r<Y+  U#oaArllgxx  '  gygxy  +  ( {Y  +  U/2)*oxxgx  =  0 

(146) 

I(Y+  U*0x9x“l[*xy  -  VYY  +  <(Y  +  1)/2)4,oxygx  "  °‘ 
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The  direction  of  the  characteristics  is  given  by 


oy/Dx  -  gy/(i-(y  +  U<t>oxgx) 


Notice  that  Eq.  (144)  cam  serve  to  express  g  by  g  . 

y  * 

It  is  convenient  to  introduce  the  arc  length  along  the 
characteristics  (to  be  denoted  by  s )  and  use  it  as  independent 
variable,  although  x  and  y  can  be  used  for  the  same  purposes, 
least  along  part  of  the  characteristics.  Then  one  has 


Dx 

Ds 


22  „ 

Ds 


i  -  (y  +moxgx 

[  (1-(Y+1) ♦ox9x) 2  +  9yl1/2 

(<l-(r+lH'ox9x)2  +  9yJ 1/>2 


and  from  Eq.  (146) 


D*x 


Ds 


d/2)92  *oxx 

[U-(yHHox9x)2  +  9211/2 


Then 


((Yf  D/2)g^oxv 
[(i-<y+i>*ox9x>4  +■  9yl 1/2 


5a  «  a  +  o  SSL. 

Ds  yx  Ds  yy  Ds* 
Finally,  from  Eq.  (145) 


D(loqH) 

Ds 


(1/2) 


<2*  1}  xx  +<^oxx^x^  ”^yy 


(147) 


at 


(148) 


(149) 


(150) 


(151) 


(152) 


(153) 
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To  determine  the  functions  gxx  and  g ■  which  occur  in  the 

last  equation,  one  differentiates  Eqs .  (146),  the  first  one  with 
respect  to  x,  the  second  one  with  respect  to  y.  One  thus 
obtains  equations  for  d9xx/Dg  and  Dg^/Ds  in  terms  of  gx,  gy,  and 
data  pertaining  to  the  basic  field.  For  the  developments  of  this 
report  these  details  are  unnecessary. 

Equations  (148)  through  (151)  represent  a  system  of  ordinary 

differential  equations  for  x,  y,  g  ,  and  g  .  The  unknowns  g,  H, 

x  y 

and  auxiliary  quantities  can  subsequently  be  found  by  quadratures . 
For  a  flow  field  given  by  <Ji  and  a  starting  line  in  the  x,  y  plane, 
one  then  can  determine  asymptotic  solutions  in  the  form  Eq.  (143) . 

For  the  usual  formulation  of  far  field  conditions  at  subsonic 
flows,  one  assumes  that  the  perturbations  originate  from  a  point. 
There  it  can  be  justified  by  the  observation  that  at  a  sufficient 
distance  from  the  origin  outgoing  perturbations  can  be  expressed 
by  a  superposition  of  particular  solutions  derived  from  Eq.  (16) . 

In  the  present  context,  such  particular  solutions  are  not  available. 
Nevertheless,  we  make  a  similar  assumption;  namely,  that  the  line 
at  which  the  far  field  conditions  are  applied  lies  at  such  a 
distance  from  the  profile  that  the  exact  point  where  the  perturba¬ 
tions  originate  does  not  matter. 

Now  we  assume  that  the  underlying  steady  flow  field  (described 
by  4>0(x,y)  can  be  approximated  in  the  far  field  by  the  dominant  term 
in  the  development  with  respect  to  the  distance.  For  flows  with  a 
free  stream  Mach  number  one,  this  term  reflects  the  profile 
displacement;  the  unsymmetry  of  the  airfoil  (including  the  effect 
of  an  angle  of  attack)  gives  rise  to  higher  order  corrections. 

The  formulae  derived  here  can  probably  be  extended  to  include  such 
terms. 

The  dominant  term  in  4>Q(x,y)  and  of  such  corrections  have  the 
form  of  similarity  solutions .  A  survey  is  found  in  Reference  8 . 

One  has  for  two-dimensional  flows 
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(154) 


7 


(1/24)c“1/5(o2  -  60  +  48) 


(159) 


r 


The  function  g  is  now  determined  by  a  similarity  hypothesis 
(rather  than  by  integrating  Eqs.  (148)  through  (152)  .  Details 
of  the  following  computations  are  found  in  Appendix  V.  Here  we 
give  only  an  overview  of  the  essential  steps. 

Let  gQ  be  an  approximation  to  g  which  satisfies  Eq.  (144) 
in  the  highest  power  of  y.  One  finds  that  gQ  is  then  determined 
from 


(Y+  1)  <b  g2  -  2g  -  g2 
'  Yox’ox  ’ox  ’oy 


=  0. 


(160) 


This  equation  arises,  if  one  omits  the  term  $  in  Eq.  (140) . 
One  usually  calls  this  the  low  frequency  approximation.  We  set 

gQ  =  y"1  y6/5a-3/5Yo(cr)  .  (161) 


Substituting  Eq.  (161)  into  Eq.  (160) ,  one  obtains 

i 

y-2y2/5o-1/  5  ( a+  ( 4/3))”1|o  ( or-  ( 16/3) )  (dYQ/da)  2-  (dYQ/da)  a  ( 2Y-  ( 20/3 ) ) 

+Yo(”3Yo  +  4)}  =  °-  (162) 

The  expression  within  the  braces  can  be  regarded  as  a  quadratic 

equation  for  dY_/da.  One  notices  that  the  coefficient  of 
2  ° 

(dY0/da)  vanishes  for  a  =  16/3.  The  curve  c  =  16/3  represents 
in  the  xy-plane  the  so-called  limiting  characteristic  of  the  far 
field.  (Perturbations  originating  downstream  of  the  limiting 
characteristic  will  never  reach  the  sonic  line  and  therefore  have 
no  effect  on  the  subsonic  field.)  In  general,  one  of  the  roots 
of  this  quadratic  equation  will  be  infinite  at  the  limiting 
characteristic.  The  solution  which  is  smooth  throughout  the 
flow  field  is  immediately  found  by  inspection 

Y0  -  4/3  (163) 


58 


Thus  from  Eg.  (161) 


go  =  U_1(y2/a)3/5(4/3) .  (164) 

The  approximation  of  the  next  order  is  obtained  by  setting 

g  -  gQ  +  (i65) 

One  assumes  that  for  y  sufficiently  large,  g^  is  sufficiently 
small  so  that  nonlinear  terms  can  be  disregarded.  A  hypothesis 
of  a  similar  kind  can  probably  be  used  to  take  higher  order 
terms  in  the  development  of  d>  into  account.  One  obtains  from 
Eq.  (144) 

(<Y+  D+oxgox  -l)glx  -  goy  gly  +  (1/2)  -  0.  (166) 

It  follows  from  Eq.  (161)  that  gQy  has  a  factor  y1^5.  A  comparison 
of  the  second  and  third  term  in  the  last  equation  then  suggests 
that  g^  has  the  form 

gL  -  (y2/a)2/5Y1(a) .  (167) 

This  leads  to  the  equation 

a(dYj/d a)  +  (2/3)yx  *  (l/4)/a  +  (4/3).  (168) 

The  one  solution  that  is  smooth  at  the  origin  is  given  by 

*  Yx  -  (1/2)  +  (3/20) a.  (169) 

One  thus  has  the  following  approximation  for  g 

g  -  u”1(y2/<7)3/5I(4/3)  +  A(  (1/2) +  (3/20)  0)]  (170) 
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with 


A  -  y2(c/y2)1/5.  (171) 

2 

Equation  (158)  gives  a  formula  for  (u/y  )  which  applies  in  the 
vicinity  of  a  =  0,  y  =  0.  In  the  next  step  one  would  determine 
the  function  H.  This  is  shown  in  Appendix  V  and  one  finds 
that  H  amounts  only  to  a  correction  of  higher  order  in  y. 

With  an  approximation  to  g  known,  one  can  derive  far 
field  conditions  in  the  following  manner.  At  a  given  point  of 
the  flow  field  the  characteristic  direction  can  be  computed  from 
Eq.  (147) .  Then  one  obtains  from  Eq.  (143) 

D<j>/Ds  -  (D  (logH)  /Ds  -  iojDg/Ds)  <J>=  0.  (172) 

In  principle,  D(logH)/Ds  can  be  computed  from  Eq.  (153),  but 
actually  it  is  small  of  higher  order.  Equation  (172)  is  a  local 
condition  which  relates  the  derivative  in  a  certain  direction 
with  the  perturbation  potential.  Because  of  the  simplification 
introduced,  this  relation  holds  only  in  the  far  field.  The 
potential  and  its  derivative  at  the  outer  edge  of  the  computed 
field  must  satisfy  these  conditions  in  order  to  match  with  the 
far  field.  One  still  needs  formulae  in  terms  of  a  and  y,  for 
these  are  the  variables  which  are  used  in  the  computations .  Here 
only  the  final  formulae  are  written  down;  the  derivation  is 
found  in  Appendix  V. 

It  can  be  assumed  that  for  a  given  point  (x,y) ,  the 
value  of  a  is  available.  For  constant  a,  A  defined  in  Eq.  (171) , 
decreases  with  increasing  y.  Only  first  order  terms  in  A  are 
taken  into  account.  The  following  formulae  are  related  to  the 
direction  of  the  characteristics.  In  practice  one  will  probably 
use  the  third  one 


D(logy)/D(loga)  =  (4/3) (1  +  (l/8)A(a+  (4/3)) 


(173) 


D(log(-x) )/D(logo)  = ( 2/3)  [ (1- (5/4) a)  ( 1-  (a/2) ) "1  +  (l/5)A(a+  (4/3))] 

D ( log ( -x))/D( logy )=  (1/2) [ (1- (5/4) a) (l-(a/2) ) _1  +  (l/5)A(a+  (4/3))] 

[1  +  (1/8) A(a+  ( 4/3) ) ] -1 

At  x  =  0,5  =  0  and  a  -  2.  For  this  vicinity,  Eq.  (173)  is  used  in 
the  form 

D(-x)  /D(loga)  =  (1/2)  [ {  (-x)  /(l-  (a/2)  }  (l-(5/4).a)  +  (174) 

(-x) (l/5)A(a+  (4/3))] [1+  (l/8)A(a+(4/3) ) ]~L 

The  term  (-x) / (1- (a/2) ,  which  at  x  =  0  is  undetermined,  is  found 
from  Eq.  (157)  and  (158) , 

(-x)/(l  -  'a/2))  =  u(a/y2)"2/5.  (175) 

No  difficulty  arises  for  a-^0,  y  =  0,  ifx/0.  One  will,  of 
course,  evaluate  Dy/Dx  rather  than  Dx/Dy. 

The  above  equations  are  correct  only  to  the  first  order  in  A. 
Some  second  order  terms  have  been  included,  however,  so  that  no 
jumps  occur  in  the  characteristic  slope  if  one  passes  from  the 
formula  for  D ( logy) /D( logo)  to  that  for  D(logx)/D(loga) . 

From  Eqs.  (173)  one  can  find  approximate  analytic  expressions 
for  the  characteristics.  For  the  formulation  of  the  far  field  con¬ 
ditions  they  are  not  necessary;  but  they  are  useful  for  a  comparisr~ 
with  characteristics  found  by  direct  computation;  One  obtains 

y  =  c  a4/3  +  c3/5  u2(a/4)(a-  (8/3)).  (176) 

Here  c  is  a  parameter  which  is  constant  along  a  characteristic. 
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One  needs,  in  addition,  Dg/Ds  and  (for  completeness)  D(logH)/DS. 


One  has 

Dg/Dc  =  (3g/3 a)  +  ( 3g/3y) (Dy/Dc) . 

This  leads  to 

Dg/D(loga)  =  u"1 (y2/a) 3/5  {[-(4/5)  +  A(-(l/5)  +  (9/100)o] 

(177) 

+  D ( logy) /D( logo) [ (8/5)  +  A(4/5)((l/2)  +  (3/20)a)]} 

Here  D(logy)/D(loga)  is  found  from  Eq.  (173).  One  obtains  in  the 
lowest  approximation 

Dg/Do  =]i^(4/3)a  1(y2/a)3//5 

Dg/Dy  =  u'1  (y2/a)  3/5  i178> 

and  by  substituting  Eq.  (176) 

-1/4 

Dg/Dy  y  '  . 

In  the  lowest  order  approximation  of  remarkable  simplicity 

D(logH)/D(logy)  =-(1/2)  (179) 


or 


D(logH.)/Dy  ~  y'1. 


Accordingly,  the  contributions  of  D(logH)/dy  in  the  far  field  condition 
is  small  in  comparison  to  that  of  Dg/Dy. 

Equations  (173)  serve  to  determine  the  direction  of  the 
characteristics  Dx/Dy  and  Do/Ds. 
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One  has 


Ds/Do  =  ((Dx/D0)2  +  (Oy/Da) 2) 1/2 . 

This,  then,  allows  one  to  compute  for  a  given  point  (x,y)  all  data 
needed  to  evaluate  Eqs.  (172) . 
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SECTION  IX 


UNSTEADY  PERTURBATIONS  IN  A  THREE-DIMENSIONAL  FLOW  FIELD 
WITH  A  FREE  STREAM  MACH  NUMBER  ONE 

The  treatment  of  the  three-dimensional  problem  is  entirely 
analogous  to  that  of  the  plane  problem.  At  great  distances  the  flow 
field  generated  by  a  body  of  finite  dimensions  approaches  that  of  an 
axial  symmetric  body.  The  approximation  to  the  basic  field  used  here 
is  therefore  the  axisymmetric  solution.  As  in  the  plane  case  the 
characteristics  are  lines;  for  an  axisymmetric  flow  field  these  lines 
happen  to  lie  in  the  meridian  planes.  In  dealing  with  an  axisymmetric 
body,  the  contour  shown  is  the  intersection  of  the  surface  with  a 
meridian  plane.  The  characteristics  in  the  meridian  plane  determine, 
of  course,  an  axisymmetric  surface;  but  this  surface  has  no  meaning 
in  its  own  right. 

The  following  formulae  are  completely  analogous  to  those  of 
the  preceding  section.  The  coordinate  y  stands  for  the  radius  in 
cylindrical  coordinates.  The  counterpart  to  Eq.  (140)  is 

— <y+  1)  (*ox  $x>  *  $yy  +  ?y/y  -  2*xt  -  ftt  -  0.  (180) 

The  hypothesis 


$(x,y,t)  =  $(x,y)exp(ia>t)  (181) 

yields 

"(Y+  V  dc  (4>OX  ♦X)  +  4>yy  +  +  ^  ^  °*  (182) 

One  sets 


<J>  -  h(x,y)  exp(-iujg(x,y) ) 

2 

The  terms  with  power  w  in  Eqs.  (182)  give 

<y+  i)*oxg*  -  2gx  -  ^  +i-o. 


(183) 


(184) 
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An  approximation,  H(x,y)  for  H(x,y)  is  obtained  with  the  power  cu 
C<Y+U<f>ox9x  “  1]  (Hx/H)  -  gy  (Hy/H)  +  (Y+l)/2)3(<J>oxgx)/3x 
-  (l/2y) 3(ygy)/3y  =  0. 

Equation  (184)  agrees  with  Eq-.  (144) .  Eqs.  (147)  through  (152) 
therefore  can  be  taken  over  immediately.  The  equation  corresponding 
to  Eq.  (153)  assumes  the  form 


DloqH  .  . 

DS  [  (1-  ( Y+  l)<foxgx)J  +  )V2 


(185) 


Next,  the  expression  in  <(>o  which  dominates  at  a  great  distance  is 
introduced 


with 


Again , 


<p  =  a3  -2/7r(r) 


-  ^-1 

5  *  U  5 


C  *  (y+  1)  _/,3xy"4//7. 


U  *  (y+  i)1/3a 

?  *  u”1  xy-4^7 


(186) 


(187) 


35/3X  =  u  ^y“4/^7 
35/  y  -  -(4/7)y“1T. 


(188) 
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The  following  closed  form  of  f  has  been  given  by  Randall 


Z  -  c“2/7(2a-l) 

I  =  -(8/9)a1/7(-2a2  +  3a  +  6) 


(189) 


By  steps  analogous  to  the  case  of  plane  flows,  one  obtains 

g  =u”1  (y2/c) 5/7 [(3/10)  +  A  ( ( 1/2)  +  (5/7) a)]  (190) 

with 

A  -  y  2  (cr/y2)  3/7 .  (191) 

D ( logy) /D( logo)  =  (6/5) (1  +  (4/3)A(o+  (1/5)]  (192) 

D (log (-x))/D (logo)  =  (2/5) [ (l-7a) (l-2a) “1+(16/7) A(a+(l/5) ) ] 

D ( log ( -x) ) /Dlogy)  »  (1/3) [ (l-7a) (l-2a) "1+  (16/7) A(a+  (1/5))] 

[1  +  (4/3) A(a+  (1/5)] 

One  obtains  as  approximation  to  the  shape  of  the  characteristics 
y  —  c  +  y 2  c1/7  (4o  -  (8/15) 

Dg/D(loga)  'y'1 (y2/a) 5/7  {[-(3/14)  +  Ay2 (- (1/7)  +  (25/24)a)] 

+  D (logy) /D(loga) [(3/7)  +Ay2((2/7)  +  (20/49)a)]} 
One  obtains  in  a  first  approximation 
Dg/Do  -  y_1(y2/0)5/7  a"1  (3/10) 

Dg/Dy  -  y"1(y2/0)5/7y"1(l/2) 
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Finally,  in  the  lowest  order  approximation 


D(logH)/D(logy)  =  -1. 
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SECTION  X 


VALIDITY  OF  THE  LINEARIZED  APPROXIMATION  FOR  TIME 
DEPENDENT  PERTURBATIONS  AT  MACH  NUMBER  ONE 

It  is  conceivable  that  the  velocities  caused  by  the  unsteady 
perturbations  decrease  more  slowly  as  one  goes  to  infinity  than 
the  deviation  of  the  velocities  in  the  basic  flow  from  the  free  stream 
velocity.  These  velocities  and  the  pertinent  sound  velocity  determine 
the  propagation  of  perturbations.  We  have  assumed  that  the  unsteady 
perturbations  are  small  and  that  their  effect  on  the  speed  of 
propagation  of  perturbations  can  be  disregarded.  But,  if  the 
unsteady  contributions  do  not  stay  small  in  comparison  to  those  of 
the  steady  flow,  then  the  present  analysis  is  not  valid  at  great 
distances,  although  it  might  be  valid  in  an  intermediate  range. 

We  carry  out  this  discussion  first  in  the  two-dimensional  flow. 
There  one  finds  from  the  similarity  solution  (Eq.  (154))  together  with 
the  second  of  Eqs.  (157)  that  along  a  line  a  =  const 

♦xo^  y“2/5  .  (193> 

For  the  superimposed  unsteady  flow  one  obtains  as  dominant  terms, 
respectively,  in  and  £ 

=  wH 

and 


(194) 


This  shows  that  in  a  perturbed  flow,  tp  dominates  4>  at  a  large 

distance.  One  notices  that  for  the  basic  flow,  one  has  the  power 

y_2//^  in  expression  <J>x,  and  y  +2//^  in  the  expression  4>x  for  the 

superimposed  perturbations.  This  makes  the  investigations  carried 

out  in  this  section  necessary.  The  question  whether  at  infinity 

the  perturbations  can  still  be  regarded  as  small  in  comparison  to 

the  basic  field  cannot  be  dismissed.  Decisive  is  the  behavior  of  the 

-1/2 

function  H.  Along  a  characteristic,  H  is  given  by  ft  or  y  ,  but 
characteristics  differ  from  lines  a  =  const.  According  to  Eq.  (171) 
one  has  along  a  characteristic 

-4/3 

ya  -v  const. 


Accordingly , 


H  =  y"1/2f(ya"4/3)  (195) 

where  f  is  used  as  a  generic  expression.  It  has  nothing  to  do  with 
the  functions  f  used  to  describe  the  basic  field.  Some  information 
about  f  is  obtained  by  considering  the  function  H  in  the  vicinity 
of  the  negative  x  axis.  It  is  reasonable  to  assume  that  in  the  x,y 
plane,  H  is  a  smooth  function.  If  one  approaches  the  negative  axis 
along  a  line  (-x)  ■  const,  then  by  Eq.  (175) 

-2 

ay  *  const, 

but  H  is  likely  to  behave  in  the  vicinity  of  the  negative  x  axis  as 
a  power  of  (-x) . 

Expressing  the  assumption  that  on  this  path  H  remains  smooth, 
one  obtains  from  Eq.  (182) 

H  -  y“1/2f(ya"3/4H  (ay"2)6  (196) 

Obviously,  the  function  f  must  have  a  singularity  at  the  x  axis. 

We  write, 
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-3/4.  ,  -3/4, ct~,  -3/4, 

f  ( ay  )  =  ( ay  )  f  ( ay  )  . 


If  one  moves  to  infinity  along  a  line  a  =  const,  then  a y”' 4/1/3  -►  0. 

In  other  words,  in  this  limiting  process  the  argument  of  f  goes  to 
zero,  throughout  the  flow  field,  one  can  choose  a  in  such  a  manner 
that  f  then  assumes  a  finite  limiting  value.  Thus,  one  obtains  from 
Eq.  (196) 


-1/2,  -3/4,  A  .  -1/2, 6 

y  (ay  )  ^  (ay  )  . 


The  powers  of  a  and  y  on  the  left  and  right  must  match 

^  a  =  -26 


a  =  6 


Hence 


a  =  2/5 


and 


H  »  (a/y2) 2/5  f(ay~3/4) 


One  finds  accordingly  that  for  superimposed  perturbations 


I  =  _2/5  -4/5  ^  -2/5 

vx  y  90X  y 


Accordingly,  the  perturbations  have  the  same  dependence  upon  y  along 
a  line  a  =  const  as  the  basic  flow. 

The  same  result  is  found  for  the  three-dimensional  flow.  One 
has  in  the  basic  steady  flow  along  a  line  a  *  const 


*x 

*x 


-V  y-6/7 

^  y+6/? 
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H  =  y  Iffy/o^5)  =  const. 


For  lines  -x  =  const.  One  has  again  (a/y-2)  =  const. 

Then  one  has  in  a  corresponding  computation 

a  =  8  -  (6/7) 

and  H  ^  y”^2/7 

along  a  line  a  =  const. 

Consequently,  the  basic  flow  and  the  superimposed  unsteady 
perturbation  change  with  the  same  power  of  y  along  a  line  a  =  const. 
If  the  unsteady  perturbations  are  small  at  an  intermediate  distance, 
they  will  be  small  also  at  large  distances. 


SECTION  XI 

COMPARISON  OF  EXACT  AND  ASYMPTOTIC  EXPRESSIONS  FOR  g 

The  analytical  formulae  for  g  derived  in  Sections  VIII  and  IX 
are  valid  only  at  a  sufficient  distance  from  the  origin.  We  show 
here  a  comparison  of  these  analytical  results  with  data  obtained  by 
a  numerical  integration  of  the  characteristic  equations.  Of  course, 
these  computations  need  not  be  performed  if  one  merely  wants  to 
apply  the  far  field  conditions . 

For  the  plane  flow  with  free  stream  Mach  number  one,  the 
equations  to  be  integrated  are  Eqs.  (148)  to  (152).  In  these 
equations  the  underlying  steady  flow  enters  via  the  function  <J>q. 

Here  the  asymptotic  expressions  for  the  steady  flow  are  used.  One 
has  according  to  Eqs.  (A. 50) 

(Y+  1)  <)>ox  =  y2(a/y2)1/5((l/4)a  -  1)  (204) 

Hence  with  Eqs.  (A. 49) 

(Y+  1Hoxx  =  ^/y2)  3/5  (a- (2/3))  (a  +  (4/3))"1  (205) 

(Y+  1>*  =  U2y  (a/y2)  6/5  (-(1/2)  a+  (4/3)  (a  +  (4/3))"1 

|i  is  a  constant  determined  by  the  underlying  steady  flow. 

In  these  expressions  the  variable  a  is  encountered.  It  is 
expressed  in  terms  of  x  and  y  (which,  for  the  integration  along 
the  characteristics,  are  dependent  variables)  by  Eqs.  (157)  and 
(158)  .  To  determine  a  for  given  x  and  y,  one  first  evaluates 

1  -  y_1xy_4/5  (206) 

and  then  determines  a  from  Eq.  ((159),  i.e., 

?  -  a"2/5((c/2)  -  1) .  (207) 
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This  can  be  done  by  a  Newton  procedure.  One  obtains  a  correction 
A a  to  an  approximation  aQ  from 


C  ~^ao) 
(d;/da)/a=a 

o 


(208) 


C(aQ)  and  d^/da  are  found  from  Eqs .  (207)  and  (A48),  respectively. 

The  method  deteriorates  in  the  vicinity  of  the  negative 
x-axis.  There  There  one  iterates  as  follows 


a(n+1)  »  (-I)“5/2(l  -  (avll,/2) 


(n)  ,,.5/2 


The  superscript  n  refers  to  the  iteration  number.  In  principle, 
one  must  prescribe  g  along  some  initial  curve.  Let  this  curve  be 
given  by 


x  =  x (y)  or  y  =  y(x) 

and  consider  a  fixed  point  (xQ,yo)  of  this  curve.  Therefore,  one 
has  for  the  system  (Eqs.  (148)  through  (152)),  initial  values  for 
x  =  xQ,  y  =  yQ,  and  g  =  g(xQ,yo)  •  One  can  also  assign  an  initial 
value  for  H.  To  find  initial  condition  for  gx  and  g^,  one  forms 
the  derivative  of  g  along  the  initial  curve 

dg/dx  *  g  +  g  (dy/dx) .  (209) 

x  y 

Here  dg/dx  and  dy/dx  are  known.  Combining  this  relation  with  the 
partial  differential  equation  for  g,  Eq.  (144)  ;  namely, 

(Y+  D<!>ox9x  "  2gx  ’  *y  +1  =  0'  (210) 

one  can  determine  g  and  g  . 

x  y 
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This  is  the  general  procedure.  In  our  computation,  we 
assumed  the  initial  curve  to  be  a  small  circle  around  a  point 
(xQ,yQ) ,  whose  radius  goes  to  zero.  From  this  circle,  characteristics 
go  out  in  all  directions.  In  this  case  one  has  the  same  value  of 
<f>  for  all  starting  points  of  the  characteristic.  Moreover, 

Eq.  (209)  can  be  disregarded.  If  one  assumes  that  g  =  const  along 
this  small  circle,  then  dg/dx  =  0.  For  any  choice  of  gx  and  gy, 
one  can  find  a  direction  (dx/dy)  along  a  circle  for  which  Eq.  (209) 
is  satisfied. 

As  initial  condition  for  the  integration  of  the  system 

(Eq.  (148)  to  (152)), one  then  has  x  =  xq,  y  =  y  ,  g  =  0.  The 

derivatives  g„  and  g  must  be  chosen  in  such  a  manner  that  Eq.  (210) 

x  y 

is  satisfied.  To  obtain  characteristics  whose  direction  is  reasonably 
spaced,  we  have  proceeded  in  the  following  manner.  Equation  (210) 
is  rewritten  in  the  form 

(-<Y+lHx<Vy0)>'1[-<Y+l>*x<xo,yo)gx+l]2+g^  =C-(y+D  *x<x0,y0) +1)  1 

[-(Y+l>*x(x0,y0) l'1 

This  equation  is  satisfied  if  one  sets 

gy  =  r  sine 

[>(Y+l)<|)x(xo,yo)gx+l]  [-(y+1)  (J)x(x0,y0)  l*1  =  r  cose 


with 


r  =  [-  (y+1)  <f> v ( x _ , y _ )  +1  )1/2]  [-(Y+l)<M*^Yjr1/2. 


Then, 

9X  "  tr(-(Y+l)<^x(x0y0)  )cose-l]  [~(Y+l)<|>x(x0,y0]"1 

The  initial  direction  of  the  characteristic  is  found  from  Eq.  (147) 
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Dy/Dx  =  (- (y+  l)^x(xo,yo)-1/2tg0 


The  axisymmetric  case  is  treated  in  an  analogous  manner .  In 

this  case  the  point  (xQ,y  )  is  chosen  on  the  negative  x  axis.  If 

one  chooses  it  away  from  the  axis  of  symmetry,  then  one  must 

assume  that  the  initial  surface  is  an  infinitely  thin  ring,  whose 

centerline  intersects  the  meridian  plane  at  the  point  (xQ,y  ) .  In 

this  case  the  characteristics  (line)  will  lie  in  the  meridian  plane. 

With  slight  modifications,  one  might  also  treat  the  case  where 

the  initial  surface  is  a  small  sphere  with  center  at  (x  ,y  ,z  =0) . 

o  o  o 

In  this  case,  one  would  obtain  a  two-dimensional  family  of 
characteristics,  which  in  general  are  not  confined  to  a  meridian 
plane . 

The  necessary  formulae  are 


(Y+l)<f>  =  U2(a/y2)  3/7  (8/3)  (<T-1) 

OX 

(Y+D^oxx  =  v  (°/y2)  5/7  ( (8/3)  a  -  (4/5))  (o+  (1/5))’1 
(Y+1)*oxy  =^2  U6/3)  y(a/y2) 10/7  (-a  +(4/5))  (a  +  (1/5)) 
1  =  u^x  y'4/7 
1  =  a“2/7(2 a  -1)  . 


-1 


The  initial  conditions  are  chosen  in  the  same  manner  as  in  the  two- 
dimensional  case. 

??he  computations  have  been  carried  out  for  y  =  1.  (See 
Eqs.  (54)  and  (186) .)  Then,  according  to  Eq.  (156)  and  (187) 


y 


(Y+  1) 


1/3 


We  have  chosen  y  -  1.4. 
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Figure  5  shows  in  a  plane  flow  with  the  free  stream  Mach 
number  one  a  set  of  computed  characteristics,  which  start  at  some 
point  of  the  negative  x-axis  and  also  the  corresponding  one-term 
approximation.  Figure  6  shows  the  corresponding  two-term  approxima¬ 
tion.  The  axisymmetric  counterparts  are  shown  in  Figures  10  and  11. 
These  curves  correspond  to  the  straight  lines  in  Figure  2.  Some 
of  them  start  in  the  upstream  direction  but  eventually  all  of  them 
are  swept  downstream.  The  curves  have  been  terminated  at  the 
limiting  characteristic  belonging  to  the  basic  flow  field.  The 
approximate  curves  have  been  chosen  in  such  a  manner  that  they 
coincide  with  the  exact  curves  at  the  outer  points.  In  practice 
exact  characteristics  are  not  available.  In  the  application  of 
the  asymptotic  far  field  conditions  the  slope  of  the  characteristics 
is  needed.  The  asymptotic  two— term  expression  gives  a  satisfactory 
approximation  to  these  slopes. 

Figures  7  (for  plane  flows)  and  12  (for  axisymmetric  flows) 
show  one-  and  two-term  approximations  for  the  wave  fronts  given 
by  lines  g  =  const.  They  are  drawn  for  g  =  10,  20,  and  30.  One 
sees  that  the  additional  term  in  the  approximation  makes  an 
essential  difference.  Figures  8  (for  plane  flows)  and  13  (for 
axisymmetric  flows)  show  computed  characteristics  and,  marked  by 
asterisk,  computed  values  of  g  (g  =  10,  20,  30)  and  also  the 
asymptotic  two-term  approximations  for  g.  The  computed  values  are 
rather  well  represented  by  their  asymptotic  curves,  especially 
if  one  takes  into  account  that  in  the  determination  of  g  a 
constant  remains  open.  The  curves  correspond  to  the  circles  in 
Figure  2.  One  no  longer  encounters  a  piling  up  of  the  wave  fronts, 
although  they  come  closer  together  at  a  large  distance  from  the 
origin.  Figures  9  show  for  plane  flow  corresponding  curves  for 
waves  starting  at  different  points.  The  asymptotic  representation 
of  the  curves  g  =  const  remains  the  same.  One  sees  that  it  is 
not  too  important  where  the  waves  originate. 


No  counterparts  for  the  three-dimensional  case  have  been 
drawn  because  the  characteristics  for  points  off  the  negative  x 
axis  will,  in  general,  not  be  curves  that  lie  in  the  meridian 
planes. 

One  recognizes  that  the  asymptotic  expressions  developed 
above  can  indeed  be  used  to  formulate  far  field  conditions, 
provided,  of  course,  that  the  boundary  of  the  computed  field  is 
not  too  close  to  the  profile. 
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APPENDIX  I 

SOME  REMARKS  ABOUT  THE  HYPOTHESIS  EQ.  (22) 


The  hypotheses  Eqs.  (22)  and  (12)  have  a  similar  effect.  In 
each  case  a  factor  is  split  off  from  the  expression  for  the 

2  A 

potential.  The  factor  exp(io)M(l-M  }  x)  in  Eq.  (12)  generates  a 
field  with  equal  waviness  in  all  directions.  Because  of  the 
interpretation  of  g(x,y)  given  by  Eq.  (36),  the  factor  exp (-iwg (x,y) ) 
in  Eq.  (22)  anticipates  the  waviness  of  the  flow  field  completely 
but  only  for  outgoing  waves.  The  original  differential  equation 
admits  incoming  as  well  as  outgoing  waves.  This  holds  also  for 
the  differential  equation  for  h  which  one  obtains  by  setting  the 
remaining  terms  in  Eq.  (23)  equal  to  zero.  To  illustrate  the 
effect  of  the  transformation  in  Eq.  (.22)  we  consider  the  problem 
for  three-dimensional  perturbations  in  air  at  rest.  There  the 
discussion  is  particularly  simple  because  of  the  availability  of 
closed  solutions.  Accordingly,  we  consider  the  partial  differential 
equation 


4>rr  +  (2/r)iJ>r 


0 


The  hypothesis 


< p  =  <Mr)  exp(iu)t) 


then  leads  to 

<j>rr  +  (2/r)$r  +  w2$  =  0 

Particular  solutions  are  given  by 

a  —  1 

4>  ■  r  exp(±iwr)  (A.l) 

(and  also  by  derivatives  of  this  expression  with  respect  to 
x,y  or  z) 
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If  one  sets  in  analogy  to  Eq.  (22) 


$  =  h(r)  exp  (-iuig  (r) ) 
then  one  has  immediately 


(A. 2) 


glr)  =  r 

The  two  particular  solutions  for  h  corresponding  to  those  in  Eq.  (40) 
are  given  by 


h(r)  =  r"1 

(A.  3) 

and 

h(r)  =  r_1  exp(i2u>r) 

(A. 4) 

In  the  expression  for  outgoing  waves  (Eq.  (A. 3)),  the  waviness  has 
vanished  while  that  for  incoming  waves  (Eq.  (A. 4))  has  waves  which 
possess  half  of  the  original  wave  length. 

This  has  interesting  consequences  for  a  numerical  approach. 

If  one  knows  that  only  outgoing  waves  are  present,  then  the  function 
h  is  smoother  than  the  function  <p.  The  mesh  for  computing  h  can, 
therefore,  be  coarser  than  that  for  $.  But  one  ought  to  be  aware 
of  the  fact  that  in  a  coarse  mesh  incoming  waves  will  be  greatly 
distorted. 

In  cases  where  the  partial  differential  equation  can  be 
solved  by  a  product  hypotheses,  so  that  ultimately  one  is  led  to 
an  ordinary  differential  equation  in  r,  one  is  tempted  to  build 
up  the  solutions  of  this  ordinary  differential  equations  from 
particular  soltuions  which  are  obtained  by  solving  initial  value 
problems.  Roughly  speaking,  the  stiffness  of  such  ordinary 
differential  equations  is  determined  by  the  largest  eigenvalue. 

If  one  solves  the  differential  equation  by  a  predictor-corrector 
method,  then  the  stiffness  causes  an  instability  of  the  procedure 
which  is  remedied  by  a  reduction  of  the  step  width.  This  happens 
even  if  the  particular  solution  which  is  being  determined  does  not 


contain  the  particular  solution  which  causes  stiffness.  The  stiff¬ 
ness  of  the  differential  equation  for  <p  is  determined  by  the  value 
of  u),  that  of  the  differential  equation  for  h  by  2u).  Of  course, 
this  discussion  is  rather  academic,  for  in  those  cases  where  a 
product  hypotheses  can  be  used  it  is  likely  that  the  ordinary 
differential  equations  can  be  solved  in  a  closed  form. 
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APPENDIX  II 

GEOMETRIC  INTERPRETATION  OF  CERTAIN  ANALYTIC  EXPRESSIONS 


First  we  determine  the  curvature  of  a  curve  g(x,y)=const. 

At  a  point  xo,yQ  we  introduce  a  local  system  of  Cartesian  coordinates 
£,  £,  where  c  is  normal  to  the  curve  g(x,y)  =  g(xQ,yQ)  and  £ 
tangential  to  it  (Fig.  14)  . 

Developing  g  with  respect  to  £  and  c  but  retaining  only  the 
terms  of  the  lowest  order,  one  obtains 

g(x,y)  =  g(xo,yQ)  +  | grad  g|c  +  g ^  \ 

Here  it  is  assumed  that  the  £  axis  points  in  the  direction  of 
increasing  g.  The  curve  g(x,y)  =  g(xQ,yQ)  is  then  given  by 

-1  £2 

C  =  “I  grad  g |  g^ 

The  curvature  of  a  curve  y  =  f(x)  is  given  by  the  familiar  formula 

-1  "  *23/2 

R  A  *  -f  /(1+f  Z)J/Z . 

It  gives  a  positive  value  for  R  if  the  region  under  the  curve 

y  =  f(x)  is  convex  (Fig.  15).  In  the  present  case  one,  therefore, 

obtains 


R"1  =  I  grad  gr1g^ 

Returning  to  the  original  x,y  system  one  observes  that 

(x  -  xQ)  *  £  cos ( £ , x)  +  ;  cos (C ,x) 

(y  -  yQ)  *  5  cos ( £ ,y)  +  C  cos ( t  ,y) 

where  (£,x),  etc.  stands  for  the  angle  between  the  respective  axes 
(here  the  x  and  the  £  axes) ,  Since  the  £  and  c  axes  are 
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perpendicular  to  each,  other  one  has 

cos(£,x)cos(£,x)  +  cos(5,y)cos(£ ,  y)  -  0 


One  then  has 


*  9xxc°s2(£,x)  +  2gxycos(C,x)cos(?,y)  +  g  cos*(£,y) 
Let,  (Fig.  16) 

(S,'x)  =  a 

•  U,y)  =  (V2  -  a) 

(C,xl  =  (ir/2  -  a) 

(£,y)  =  Or  -  ct) 


(The  signs  of  these  angles  do  not  matter  because  one  encounters  in 
the  above  formulae  only  the  cosines  of  these  angles . )  One  obtains 

2  2 
*  ^XX3in  a  “  29XySina  COS“  +  9yyCOS  “ 


and 


-1  2  222  1/2 
R  =  (g^sin  a  -  2gxysina  cosa  +  g^cos  ct)/(gx  +  gy)  x  (A. 5) 

Notice  that  in  this  formula  specific  properties  of  the  g-field 
are  not  taken  into  account.  This  is  done  presently. 

One  has 

gx  *  | grad  g | cos  a 
gy  »  | grad  g|sin  a 

Substituting  these  expressions  into  Eg.  (24)  one  finds 

| grad  g|  ■  Cl  +  M  cos  a)“^  (A. 6) 


83 


Therefore, 


gx  =  cosa(l  +  M  cosa)  ^ 
gy  =  sina(l  +  M  cosa)-1 

The  direction  of  the  characteristic  is  given  in  Eq.  (29) . 
obtains 


(A. 7) 


One  now 


Dy/Dx  =  sina/(M  +  cosa)  (A.  8) 

It  was  found  in  Section  II  that  g  ,g  and  therefore  also  a  and 

x  y 

Dy/Dx  are  constant  along  a  characteristic.  Introducing  a  parameter 
p  which,  at  the  point  under  consideration  assumes  the  value  pQ, 
one  obtains  the  following  parametric  representation  for  a 
characteristic 


x  =  x„  +  (p  -  p  ) (M  +  cos  a) 
o  o 

y  *  yQ  +  (p  “  PQ>sin  a 
Along  such  a  line  one  has  the  relation 


(A. 9) 


Dg/Dp  »  g  (M  +  cos  a)  +  g  sin  a  *  1  (A. 10) 

x  y 

Equations  (A. 9)  and  (A. 10)  lead  to  the  following  geometric 
interpretation.  We  consider  a  point  (xo,yQ)  of  the  curve  g  -  gQ. 
According  to  Eq.  (A. 10) ,  the  value  of  p,  to  be  denoted  by  p^,  for 
a  curve  g  *  g^  is  given  by 

Pi  -  p0  =  •  v 

One  has  for  the  point  x^,y^,  originating  from  (xQ,po) 

<xl-xo>  *  (P1“P0)  cosa  +  (Pj-P^M 

<y0-y0)  -  (P1“PQ) sina. 


f* 


Accordingly,  the  line  g  »  g^  arises  from  the  line  first  g  =  gQ, 

by  proceeding  in  the  direction  of  the  normal  by  a  distance 

(p^-pQ)  =  (g^-gQ)  an<*  afterwards  by  translating  .he  curve  so 

obtained  in  the  x-direction  by  a  distance  (p^-pQ)M.  The  first 

step  amounts  to  the  determination  of  a  curve  equidistant  from 

g  *  gQ  by  the  distance  (p^-p^ •  By  this  process ,  the  radius  of 

curvature  is  changed  from  R  (at(x  .y,)  to  R,  =  +  p  -  p  (at 

o  o  o  x  o  o 

the  point  corresponding  to  (x^,y^)]U  The  subsequent  translation 
leaves  the  radius  of  curvature  unchanged.  Identifying  po  with  Rq 
one  finds  that  at  the  points  of  subsequent  curves  g  =  const,  the 
parameter  p  is  identical  with  the  radius  of  curvature. 

The  analytical  derivation  of  this  result  is  somewhat  cumber¬ 
some.  One  can  proceed  as  follows.  The  fact  that  g  and  g  are 

x  y 

constant  along  characteristics  is  now  used  to  expresss  g  ,g  , 

_  i  3cx  y  y 

and  g^  in  terms  of  R  .  It  follows  from  Eq.  (A. 9)  that 


Dgx/Dp  =  gxx(M  +  cos  a) +  gxy  sin  a-  0 
Dgy/DP  *  gxy(M  +  COS  a) +  gyy  sin  a  *  o. 
Moreover,  from  Eq.  (A. 5)  in  conjunction  with  Eq.  (A. 6) 


(A. 11) 


2  2  -1  —1 
gv„  sin  a  -  2g  sina  cosa  +  g„„cos  a«  R  (1+Mcosa)  (A. 12) 

a*  xy  yy 

Equations  (A. 11)  and  (A. 12)  form  a  system  of  linear  equ'tions  for 
gxx'gxy'  and  gyy *  0ne  verifies  that 

g„„  *  r""1  sin2a/(l  +  M  cos  a)  3 

g  *  -R-1  sin  a (M  +  cosa)/(l+M  cosa) 3  (A. 13) 

xy 

g  ■  R~1(M  +  cosa)2/(l  +  M  cosa)3. 

yy 
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Hence 


g  (1-M2)  +  g  =  R-1/ (1  +  M  cos  a)  (A. 14) 

xx  yy 

In  these  expressions  R  is  the  local  radius  of  curvature. 

So  far  the  relation  between  R  and  p  has  not  been  established. 

For  this  purpose  we  express  the  derivative  of  R  ^  along  a 
characteristic  from  Eq.  (A. 14) ,  keeping  in  mind,  that  a  =  const 
along  a  characteristic 

-R_2(dR/dA)  =  (1  +  M  cosa) [(1-M2) (dgxx/D£)  +  (Dg  /DA) ]  (A. 15) 

The  second  derivatives  of  g  encountered  here  are  subject  to 
compatibility  conditions  (as  are  all  higher  derivatives) ,  which 
are  derived  from  Eq.  (28)  .  One  obtains  by  forming  derivatives  with 
respect  to  x  and  y 

[(1-M  )gx  +  M)gxxx  +  9y9XXy  + 

[U-M2)gx  +  HJ9yyx  +  gygyyy  + 

[ (1-M2) g  +  MIg  +  g  g  + 

x  ^xyx  y  xyy 

It  follows  from  Eqs.  (A. 7)  that 

[1-M2)gx  +  M]  =  (M  +  coscx)  /  (1  +  M  cosa)  . 

Then  with  Eqs.  (A. 7)  and  (A. 9) 

Ogxx/DP  +  (1  +  M  cos  a)  I(l-M2)g2x  +  g^]*  0 

Dgyy/Dp  +  (1  +  M  cos  a)  t(l-M2)g2y  +  g2]  =  0  (A. 16) 

Dgxy/Dp  +  +  M  cos  $  [  ( !-M2 )  gxx  +  SyylSxy  = 


(1-M2) g2  +  g2  =0 
'  *xx  ^xy 

a-M2)g2y  +  gyy  -  o 

(  1  M  )  gXy9XX  +  gyy9xy  “  0  . 
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0. 


This  is  a  system  of  three  nonlinear  ordinary  differential 

equations  for  g  ,  g  ,  and  g  .  It  holds  along  a  characteristic, 
xx  yy  xy 

The  previous  result,  that  p  is  identical  with  R,  which  in  the 
present  context  is  to  be  regarded  as  a  conjecture,  suggests  that 
one  particular  solution  is  given  by  Eq.  (A. 13)  with  R  replaced 
by  p.  If  this  is  correct,  then  one  has 

-2  2  3 

Dg^/Dp  =  -p  sin  a/  (1  +  M  coset) 

—  2  ^ 

Dgxy/Dp  =  p  sina(M  +  cos  a)/(l  +  M  cosa)  (A. 17) 

Dgyy/Dp  =  “p_2^m  +  cos  +  M  cos 

Equations  (A. 13)  (with  R  replaced  by  p)  are  indeed  particular 
solutions.  One  obtains,  for  instance,  by  substituting  into  the  first 
of  Eqs.  (A. 16) 

-p~2sin2a (1+Mcosa) -3+p~2 (1+Mcosa) [ (1-M2) sin4a+  sin2a(M+cosa) 2]  =  0. 

Equations  (A. 13)  give  only  a  particular  solution,  for  they  do  not 
contain  three  constants  of  integration.  Actually,  one  does  not 
need  the  general  solution,  for  at  the  initial  point  (p  =  pQ) , 

Eqs.  (A. 13)  (with  R  replaced  by  p)  assume,  of  course,  the  values 
of  gxx'gxy'  and  gyy  expressed  in  terms  of  RQ. 

So  far  we  have  identified  R  with  p  only  at  the  point  (xo,yQ) . 

The  identification  at  other  points  of  the  characteristic  is  still 
a  conjecture.  Substituting  Dg^^Dp  and  Dg^^/Dp  (Eq.  A.  16)  into 
Eq.  (A. 15)  and  remembering  that  according  to  Eq.  (53)  dp/d£=const 
along  a  characteristic,  one  obtains 

-R-2 (dR/dp) »- (1+Mcosa) p“2 [ (1-M2) sin2a  +  (M  +  cosa) 2]/ (1+Mcosa) 2 . 

-R~2 (dR/dp)  «  -p“2. 


Hence 


+  const. 
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The  constant  vanishes,  because  one  can  choose  at  the  point  (xQ,yo) 

ro  o 

This  fits  with  the  above  conjecture.  The  identification  of  p  with 
the  radius  of  curvature  is  therefore  justified. 

To  determine  H^  we  rewrite  the  first  of  Eqs.  (34)  using 
Eqs.  (53) 

D(logH1)/Dp=-(l/2) [l-M2)gxx+gyy] (gy  +( (1-M2) gx+M) 2]~1/2 [M2+2Mcosa+l] 1/2 
Hence,  with  Eq.  (50) 

DUogH^/Dp  =  -(1/2)  [l"M2)gxx  +  gyy]  (1  +  M  cosa) 
and  with  Eqs.  (A. 13)  with  R  replaced  by  p 

DdogH^/Dp  =  -  ( 1/2)  p-1 . 


Hence, 


H^  =  const  p  ^2  (A. 18) 

The  form  of  the  function  (introduced  in  Eq.  (25))  is 
suggested  by  the  asymptotic  development  of  the  Hankel  functions, 

(Eq.  (19)).  In  the  present  setting  K2  is  found  by  integrating  the 
second  of  Eqs.  (34).  With  p  as  independent  variable,  Eq.  (34) 
assumes  the  form 

DH2/dp  +  (l/2)p“1H2  +  (i/2)  Kl-M2)^^  +  H^yy]  (1+Mcosoe)  =  0  (A. 19) 

Crucial  is  the  inhomogeneous  term.  The  constant  in  Eq.  (A. 18)  may 
differ  from  characteristic  to  characteristic.  Accordingly,  we 
write 

Hx  *  f^sjp'1/2  (A. 20) 
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where  s  is  the  arc  length  of  the  initial  curve  (see  Eqs.  (47)  and 
(48)).  The  value  of  s  is  then  attached  to  the  characteristics  that 
start  at  the  initial  curve.  To  estimate  the  inhomogeneous  terms, 
we  proceed  as  follows.  Introducing  s  into  Eqs.  (A. 8)  we  write 


x  =  x(p,s)  =  xQ(s)  +  (p-pQ(s))(M  +  cosa(s)) 
y  =  y(p,s)  =  yQ(s)  +  (p-pQ(s))  sin  a(s)  . 
Inverting  this  transformation,  one  obtains 


(A. 21) 


p  =  P(x,y) 
s  =  s  ( x ,  y )  . 

One  has  the  following  relations 

3p/3x  =  (3y/3s)/D 

3s/ 3x  =  -(3y/3p)/D 

(A. 22) 

3p/3y  =  -(3x/3s)/D 

3s/3y  *  (3x/3p)/D 

with 


D  =  (3x/3p) (3y/3s)  -  ( 3y/3p) / ( 3x/3s)  (A. 23) 

Differentiating  Eqs.  (A. 21)  one  obtains 

3x/3p  =  M  +  cosa(s) 

3y/3P  *  sina(s) 

3x/3s  ■  -(3p  /3s)(M  +  cosa)  +  ( 3x  /3  ) - (p-p  (s) ) sina (3a/3s) 

o  o  s  o 

3y/ 3g  »  -(3pQ/3s) sina  +  (5yQ/3s)  +  (p  -pQ (s) ) cosa ( 3a/3s) . 
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The  parameter  s  is  the  arc  length  along  the  initial  curve,  a  the 
angle  of  its  normal  with  the  x  axis,  and  pQ  the  local  radius  of 
curvature.  Therefore 

3a/ 3s  =  p"1 
3xq/3s  =  -  sin  a 
3yQ/3s  *  cos  a 

and 

3x/3s  =  -(3pQ/3s)  (M  +  cosa)  -  (p/pQ)sin  a 
3y/3s  =  - { 3pQ/3s) sin  a +  (p/pQ)cos  a 
and  from  Eq.  (A. 22) 

D  =  (p/pQ)  cos  a  +  1)  . 

Then,  from  Eqs .  (A. 22) 

3p/3x  =  [- (pQ/p) ( 3pQ/3s) sin  a  +  cos  a] (M  cos  a +  1)”^ 

3s/3x  =  -  (pQ/p)  sin  a  (m  cos  a  +  1)“  (A. 

3p/3y  =  [“(P0/p)  (3pQ/3s)  (M  +  cosa)  +  sina]  (M  cos  a  +  l)  ^ 
3s/3y  =  (pQ/p) (M  +  cos  a) (M  cos  a  +  1) 

From  these  equations  one  finds  the  order  of  magnitude  of 
the  expressions  for  large  values  of  p.  3pQ/3s  and  pQ  are 
considered  as  quantities  of  order  1.  One  finds 


3p/3x  *  0(1);  3s/3x  =  0<p_1) ;  3p/3y  =  0(l);3s/3y  *  0(p”^) . 

The  orders  of  magnitude  are  needed  also  for  derivatives  of  the  next 
order .  One  has 


and  similar  formulae  for  the  other  derivatives. 

In  forming  derivatives  of  (3p/3x) ,  (3p/3y),  (3s/3x) ,  and 
(3s/3y)  with  respect  to  s,  the  order  of  magnitude  of  the  terms 
remains  unchanged;  in  forming  derivatives  with  respect  to  p,  the 
order  is  lowered  by  a  factor  p*"1 


32p/3x2  =  0(p_1)0(l) 

+  OllJOfp'1)  = 

0(p_1 

32p/3y2  =  0(p“1)0(l) 

+  0  (1)  0  (p"*1)  = 

0(p_1 

32s/3x2  «  0(p"2)0(l) 

+  0(p"1)0(p*1) 

*  0(p 

32s/3y2  -  0 (p""2) 0(1) 

+  0(p"1)0(p"1) 

=  0(p 

Now  we  can  estimate  the  second  derivatives  of  .  One  has 

dHj/dx  =»  (3H1/3s)  (3s/3x)  +  (3H1/9p(  3p/3x) 

a^j/ax2  =  O^/a2)  Os/3x)2+2(32H1/9s3p)  (3P/3X)  (3s/3x)  + 

+  (32H1/3p2)  (3P/3X)  2  +  OHj/as)  (32s/3x2)  +  ( 3Hj/3p)  (32p/3x2) 
2  2 

and  similarly  for  3  H^/Sy  . 


Now  according  to  Eg.  (A. 20) 


Then, 


a^/ax2  =  0(p“1/2)0(p"2)+0(p"3/2)0(i)0(p"1)+0(p 
+  0{p'1/2)0(p"2)+0(p"3/2)0(p’1) 

a^/ax2  =  o(P-5/2)  . 

Similarly, 

a^/ay2  =  0(p"5/2)  . 

-5/2 

The  inhomogeneous  terms  in  Eq.  (A.  19)  is  therefore  0(p  ) .  A 

-3/2 

particular  solution  of  this  equation  therefore  is  0(p  ) .  The 

solution  for  the  homogeneous  part  has  0(p“1//2)  ,  but  this 
contribution  can  be  incorporated  in  .  This  is  the  expected 
result  for  H2, 

To  arrive  at  an  equation  analogous  to  Eq.  (A. 14),  for  the 
three-dimensional  problem,  one  can  proceed  as  follows. 

First  one  derives  a  formula  for  the  average  curvature 

at  a  point  x  ,y  ,z  of  a  surface  g  »  const  in  the  x,y,z  space. 
c  o  o  o 

We  introduce  a  local  system  of  Cartesian  coordinates  £,  n ,  £ .  The 
5  axis  has  the  direction  of  the  normal  to  the  surface  g  *  const 
through  the  point  (x0'Y0'z0> •  The  5  and  n  axis  then  lie  in  the 
tangential  plane.  The  directional  cosines  with  respect  to  the 
x,y,z  system  are  denoted  by  cos(5,x),  cos(Cry)*  etc.  One  has 

5  2  2 
cos  ( £,x)  +  cos  (n»x)  +  cos  (cx)  ■  1. 


-5/2 


)0(1) 


Hence 

9  2  2 

cos  ( 5,x)  +  cos  (r),x)  »  sin  (z,x) 

cos2(C,y)  +  cos2(n,y)  -  sin2(c,y)  (A. 25) 

co82(C,z)  +  cos2(n#*)  *  sin2(c,z) 
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and  because  the  x,y  and  z  axis  are  perpendicular  to  each  other 


cos(  £,x)cos(  £,y)  +  cos(ri,x)cos(n,y)  =  -cos(s,x)cos(C,y) 

(A. 26) 

cos(£,x)cos(£,z)  +  cos(n,x)cos(n»z)  =  -cos  (?  ,x)  cos  (s ,  z) 
cos(  £,y)cos(£,z)  +  cos(n#y)  cos  (n/Z)  =  -cos ( c ,y) cos U , z) . 


The  coordinate  systems  are  related  by  the  equation 

x-xQ  =  £cos(£,x)  +  n  cos(n»x)  +  x,  cos(c,x 

y-yQ  =  £cos(£,y)  +  n  cos(n»y)  +  C  cos(c,y) 

z-zQ  =  £cos(£,z)  +  n  cos(n»z)  +  S  cos(?,z)  . 

Developing  g  in  the  vicinity  of  the  point  ( xQ,yo,zo ) 
of  second  order,  one  obtains 


(A. 27) 


up  to  terms 


g(x,y,z)  =  g(xo,yo,zQ)  +  |grad  g|£ 

+  g^U2/2)  +  gnn  (Ti2/2)+g^(52/2) 

+  g£n^  +  g5Sfe  +  gncn5’ 

The  intersection  of  the  surface  g  *  g(x0,y0,z0)  with  the  plane 
n  «  0*  is  then  given  by 

| grad  g|  £  +  g^C^/2)  +  g^(c2/2)  +  g^©;  -  0. 

Let  R^  be  the  radius  of  curvature  of  this  curve  at  the  origin. 
One  finds 


R"1  -  g.£5/|grad  g|  . 

Here  it  is  assumed  that  the  region  g  <  g(xQ,yo,z0)  is  convex 
and  the  center  of  curvature  lies  on  the  inner  normal  to  this 
region . 


In  a  corresponding  manner,  one  obtains  for  the  radius  of  curvature 
of  the  intersection  of  the  plane  £  »  0  with  the  surface 
g  -  g(x0,y0,z0) 

R'1  *  9nn  /lgrid  «l 


with  corresponding  rules  for  the  choice  of  the  sign. 
The  average  curvature  is  then  given  by 


'average 


(1/2) (r“X  +  Rj1)  =  (1/2) (g  +  g  ) / | grad  g| . 


(A. 28) 


We  shall  see  that  Caverage  does  not  depend  upon  the  orientation 
of  the  £  and  n  axes  (as  long  as  they  lie  in  the  tangential  plane 


We  express  (1/2)  (g^  +  g  ) 


and  are  perpendicular  to  each  other) 
in  the  x,y,z  system.  Using  the  transformation  formulae  (Eq.  (A. 27) 
one  obtains 


(1/2)  (g?c  +gnn)  -  (l/2)gxx(cos2(£,x)  +  cos2(n  ,x) )  +(1/2)  gyy  (cos2  (  £,y) 

+cos2(n,y))  +  (l/2)g  (cos2(£rz)  +  cos2(n,z)) 

zz 

+  gw  ( cos  ( £ ,  x)  cos  ( £ ,  y )  +  cos  (n,x)  cos  (n  ,y) 

+  gxz(cos(£,x)cos(£,z)  +  costn,x)cos(n,z) 

+  gy2Ccos{£,y)c0s(£,z)  +  cos(n,y) cos(n ,z) )  . 

Hence,  with  the  relations  CEqs.  CA.25)  and  (A. 26) 

(1/2)  gC5+gnn) » (1/2)  g^sin2  (?  ,x)  +U/2)  g^sin2  (c  ,y)  +  (1/2)  gzzsin2  ( £ ,  z) 

cosC?,x)cos(?,y)  -  g  cos (; ,x) cos U , z) -g„  cosU,y)cos(;,z)  . 

(A. 29) 

This  is  the  analogon  to  Eq.  (A. 14)  for  the  two-dimensional  case. 

The  directional  angles  for  the  £  and  n  axes  are  no  longer  present. 
This  shows  that  this  expression  does  not  depend  upon  the 
orientation  of  these  axes  in  the  tangential  plane. 


Again,  gx,g^,  and  gz  are  constant,  along  the  characteristics 
(see  Section  VIII)  whieh  are  given  in  parametric  form  in  &g.  (115) 
The  ;  direction  is  the  direction  of  the  normal  to  a  surface 
g  -  const.  One  then  obtains 


e 


gxx(M+cos(5,x) )  +  g  coa(c,y) 
gw(M+cosU,x) )  +  g^cosUjy) 

xy  yy 

9X2(M+C°8(i;'x)  +  9yZc©s  ( ?  ,y) 


+  gxzcos(i:,2) 
+  gyzcos(?,2) 
+  gzzcos(C,2) 


0 

0 

0 


(-M+COS (C  ,x) 
cos(c,y)  (A. 30) 
cos(C,z) . 


Multiplying  the  equations  with  the  factors  shown  on  the  right  and 
adding  the  results  to  Eq.  (A. 29)  multiplied  by  2,  one  obtains 


9«  +  9nrT  9*X(1‘M  >  +  gyy  +  9»- 


The  average  curvature  (Eq.  CA.28))  is  therefore  given  by 


(A. 31) 


Caverage  *  (1/2)  (R^  +  Rj1)  -  (1/2)  Cg^d-M2)  +  gyy  +  gzz)/|grad  g| 


In  Eq.  (115)  the  characteristics  have  been  represented  in 

terms  of  a  parameter  p.  To  determine  the  dependence  of  the 

average  curvature  upon  p,  we  use  the  geometric  interpretation  of 

Eq.  (115)  already  applied  in  this  Appendix  for  the  two-dimensional 

case.  Constructing  a  point  of  the  surface  g  ■  g^  *  const  from 

point  (xrt,y_,z  )  of  a  surface  g  ■  g  ■  const,  one  first  proceeds 
I  o  o  o  o 

in  the  direction  of  the  normal  to  the  original  surface  by  a 
distance  p1-pQ  *  g^-gQ  (because  of  Eq.  (117) ,  and  then  shifts 
.  the  surface  in  the  x-direction  by  M(p^-p0) .  In  the  first  step 
one  constructs  a  surface  equidistant  to  g  -  gQ  at  a  distance 
|  t  9l~9o'  in  ***•  *®cond  t*1®  surface  as  a  whole  is  translated  in 
the  x-direction. 

The  radii  of  curvature  R^  and  Rj  refer  to  the  lines  of 
intersection  of  a  surface  g  *  const  with  two  planes  l  ■  const 
and  n  -  const,  which  are  perpendicular  to  each  other  and  contain 

i 
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the  normal  to  this  surface.  Among  the  possible  orientations  of 

these  planes ,  there  is  one  where  and  have  their  extrema. 

One  then  ohtains  the  principal  radii  of  curvature.  For  these 

orientations  of  the  planes  the  normal  to  the  surface  g  *  const 

remains  within  the  planes  even  at  points  of  the  surface  which 

deviate  from  (x  ,y  )  by  a  first-order  distance.  One  then  obtains 
o  o  o 

points  of  the  surface  g  *  g^  by  proceeding  from  a  point  of  the 
line  of  intersection  of  g  *  gQ  with  5  =  0,  say,  to  the  "corresponding" 
point  of  the  surface  g  =  g^  by  proceeding  within  the  surface 
5=0.  In  determining  the  radius  of  curvature  of  the  intersection 
of  the  surface  g  =  g^^  with  the  plane  5  =  0,  we  can  confine  our 
attention,  to  this  plane.  If  the  radius  of  curvature  of  the  line 
of  intersection  of  the  surface  g  =  gQ  is  given  by  *  Rio'  one 
finds  for  the  corresponding  radius  of  curvature  R^,  at  q  =  g^ 

R11  "  R10  +  (pl  “  po)* 

If  one  sets  pQ  =  0,  then  one  obtains  for  the  average  curvature 
(Eq.  (A. 31)) 

Average,  1  ’  +  f'1  +  (R20  +  R  >'X> 

""average,  1  =  {  U/2)  I  Cl-M2),^  +  gyy  +  g„]/|grad  gl}, 

=  (1/2)[(R10  +  p)"1  +  (R20  +  p)-1] 

where  R1Q  and  R2q  are  the  principal  radii  of  curvature  at  the 
point  x  .y  .z  .  Only  for  large  values  of  p  will  this  expression 
behave  as  p”  . 
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APPENDIX  III 


VERIFICATION  OF  EQUATION  (67) 

It  is  practical  to  rewrite  Eqs.  (66)  and  the  next  equation. 

Let 


Furthermore 


Rl  R2  "  r3  = 

2 

=  2 -  u>  <j> . 


X  =  (1/2)  (♦„  ♦  *yy> 

Y  -  U/2)  <♦„  -  4.yy> 

Z  -  ^xy * 

The  following  trigonometric  identities  are  used 

cos2B  =  (1/2) (1  +  cos  28) 
cos2y  ■  (1/2)  (1  +  cos  2y) 
sin8sinY=  (1/2) (cos(B-Y)  -cos(8+y)>. 


(A. 32) 


(A. 33) 


(A. 34) 


For  later  use  we  note  some  further  identities 

cos28+cos2y  =*  l+(l/2)  (cos  20+  cos  2y)  “l+cos (B+y) cos (8-y) 

cos28+cos2y  »  2cos(8+y) cos(B-y) 

2  2 

cos  8-cos  y  ■  (1/2) (cos  28-cos  2y)  =  -sin(8+Y) sin(8”Y)  (A. 35) 
cos2b«cos2y  -  -2sin(B+Y) sin(B-Y) 

»in2S-sin2Y  *  2sin(8-Y) cos(8+Y) 
sin2g+sin2Y  *  2sin(0+Y) cos (8-Y) 
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Equations  (66)  then  assume  the  following  form 


X  +  Y  cos20  +  Z  sin20 
X  +  Y  cos2y  +  Z  sin2Y 
X  cos(B-y)  +  Y  cos(8+Y) 
X(2-M2)  -  Y  M2  =  R4 


+  Z  sin( S+Y) 


(A. 36) 


Equation  (67) ,  which  is  to  be  verified,  now  appears  in  the  form 

(1-M2sin2y)  +R2  (1-M2sin20)  +R3  [-2cos  (  0~y)  -M2  (cos  ( 0+y)  -cos  ( 0-y) )  ] 

(A. 37) 


-R4sini 


(b-y) 


0 


To  carry  out  this  verification  one  must  show  that  the  scalar  product 
of  the  vector  formed  with  the  left  hand  sides  of  Eqs.  (A. 36)  and 
the  vector  formed  with  the  coefficient  of  R^  through  R4  in  Eqs. 

(A. 37)  vanishes.  Accordingly,  one  must  evaluate  the  scalar  product 


X 

+Y  cos2B 

+Z 

sin2  0 

2  2  3 

1  -  M  sin  y 

X 

+Y  cos2y 

+Z 

sin2y 

1  -  M2sin20 

Xcos (0-y) 

+Y  cos (0+y) 

+Z 

sin (0+y) 

•  . 

{-2  cos (B-y)  j 

-M2 (cos ( 0+y) -cos ( B-y) ) 

X(2-M2) 

-YM2 

-sin2 (B-Y) 

2  .  .  2 
This  expression  is  linear  in  M  .  The  terms  not  containing  M 

give 

X[1+1-2cos2(8-y) -2sin2(6-Y) +Y[cos2g  +  cos2y-2cos (B+y) cos (B-y) ) 


+Z[sin20  +  sin2y  -  2sin(B+Y)  cos  (0-y)  1  • 

The  coefficients  of  Y  and  Z  vanish  because  of  Eqs.  (A. 35) .  The 
vanishing  of  the  coefficient  of  X  is  self  evident. 


P" . p  11  ' ■  1  .  IIILBIU  'lu  1  1 1  ■"  ■  ■  I  ^ 

I  1 

'i 

2 

Next  we  evaluate  the  terns  with  H  .  It  suffices  if  one 
2 

sets  M  -  1.  The  second  vector  then  becomes 

2  2  2  + 

{cos  y,  cos  Q,  -  [cos (8+y)  +  cos(e-y)],  -sin  (8~y)]  . 

One  obtains  for  the  coefficient  of  X 

2  2  2  2 
cos  y+  cos  8  -  cos  (8-y) -cos (8~y) cos (8+y) -sin  (8-y) • 

It  vanishes  because  of  the  first  of  Eqs.(A.35). 

The  term  with  Y  is  given  by 

2  2*2  2 
Y {cos (28) cos  y  +  cos (2y) cos  8  “  cos  (8+y) -cos ( 8+y) cos (8-y) +sin  (8-y) 1 

=  Y  [cos (28) (1/2) (cos (2y)  +  1)  +  cos (2y) (1/2) (cos(28)+l-cos2(e+Y) 

-cos (S+Y) cos (8-Y)  +  sin2 (8-y)l. 

Next,  with  the  second  of  Eqs.  (A. 35) 

»Y  {cos (28) cos2y  -cos2 (8+y)  +  sin2 (8-Y)] 

Now  cos28  and  cos2y  are  expressed  in  terms  of  8+y  and  8-y 
»Y [cos2 (8+y) cos2 (8-y) -sin2 (8+Y) sin2 (8-y) -cos2(6+y) +sin2 (6-y) ] 

»Y [-cos2 ( 8+y) sin2 ( 6-y)  +  sin2 (8-y) cos2 (8+y) ]  *  0. 

The  term  with  Z  is  given  by 

Z[sin(28)cos2Y  +  sin(2Y)cos2B  -  sin(8+Y) (cos(6-Y)+cos(8+Y)) ] 

-Z  [sin(28) (1/2) (cos(2Y)+l)  +  sin(2Y)(l/2) (cos(28)  +  1) 

-sin ( 8+Y) cos ( 8-Y) -sin ( 8+Y) cos ( S+Y) ] . 


99 


APPENDIX  IV 
CHARACTERISTICS 


We  compile  here  the  leading  ideas  of  the  theory  of 

♦ 

characteristics  applied  to  the  present  case.  The  differential 
equation  under  consideration  is  Eq.  (5),  with  $  replaced  by  4>. 

(1-M2)  <p  +  <p  -2M<p  .  ~  <Jk.  =  0  (A. 38) 

Yxx  Yyy  yxt  ytt 

We  consider  some  surface  in  the  x,  y,  t  space  and  assume  that  it 
is  oriented  so  that  it  can  be  parametrized  by  x  and  y.  Let  this 
surface  be  given  by 


t  =  f (x,y) 


(A. 39) 


We  assume  that  <jj  =  U(x,y)  ,  <|>  =  V(x,y)  and  <j>.  =  W(x,y)  are  known 

x  y  t 

on  this  surface.  Obviously,  these  quantities  cannot  be  independent. 
If  the  surface  S  were  given  by  t  =  const,  then  in  order  for  a 
potential  to  exist,  one  must  have  Uy  =  Vx.  To  derive  a  corresponding 
relation  for  a  general  surface,  we  consider  the  following  system 
of  equations 


10  0 
0  10 
0  0  1 
0  0  1 
0  0  0 
0  0  0 


f  0  0 

4* 

u 

X 

Yxx 

x 

0  f  0 

4> 

V 

y 

yy 

y 

f  0  0 

<() 

u 

y 

Yxy 

y 

0  f  0 

— 

V 

X 

Yxt 

X 

1  0  f 

ip  , 

W 

X 

ryt 

X 

0  1  f 

<p.. 

W 

y_ 

Ytt 
—  — 

_  y_ 

(A. 40) 


The  arguments  of  the  derivatives  of  <j>  are  x,  y  and  f(x,y)  =  t.  The 
arguments  of  the  derivatives  of  f,  U,  V,  and  W  are  x  and  y. 
Premultiplying  Eq.  (A. 34)  by  the  matrix 
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The  left  side  of  the  third  and  fourth  equations  in  this  system  are 
identical.  Therefore,  the  right  sides  must  be  identical  also. 


U  -W  f  =  V  -W  f 
y  x  y  x  y  x 


(A. 41) 


In  deriving  the  characteristics  condition,  one  will  omit  either 
the  third  or  fourth  equation,  but  add  the  partial  differential 
equation  for  Eg. (A. 38),  which  so  far  has  not  bee  used.  Accordingly, 
one  considers  the  system 
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etc .  amounts 


The  systematic  elimination  of  the  unknows  <J>  » 

xx  yy 

to  premultiplication  of  the  system  by  the  following  1  by  6  matrix: 
[-(1-M2),  -1,  0,  2M+(1-M2) f  ,  f  .  1]. 

x  y 


This  yields  the  relation 


(A. 42) 


[(2M+(1-M2)fx)fx  +  f2  -13<J>tt  =  -(l-M2)Ux-Vy+(2M+(l-M2)  fx)Wx+fyW=0 


For  a  characteristic  surface,  the  derivatives  (among  them 
4>tt)  are  not  uniquely  determined.  Along  characteristic  surfaces 
discontinuities  in  the  second  (and  higher)  derivatives  may,  therfore, 
be  encountered.  Furthermore,  derivatives  within  the  surface  cannot 
be  chosen  independently.  This  can  be  recognized  in  the  last  equation. 
One  has  a  characteristic  surface,  if  the  factor  of  $  vanishes 

( 1-M2)  f 2  +  2Mfx  +  f 2  -  1  =  0  (A. 43) 


But  then  the  right  hand  side  must  also  vanish  and  one  obtains  the 
compatibility  condition 

- ( 1-M2) U  -V  +  ( 2M+ ( 1-M2) f  ) W  +  f  W  *  0.  (A. 44) 

x  y  jk  y  y 
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APPENDIX  V 


THE  TWO  DIMENSIONAL  UNSTEADY  FLOW  FIELD  AT  MACH  NUMBER  ONE 

DETAILED  COMPUTATIONS 


In  the  main  text  the  general  course  of  the  computations  is 
shown.  The  details  are  rather  cumbersome,  although  certain 
expressions  which  at  the  beginning,  have  a  rather  unwidely  form 
can  be  considerably  simplified.  The  present  appendix  gives  these 
computations  in  sufficient  detail  to  allow  a  sceptical  reader  to 
check  the  results.  The  function  g  Is  determined  from  the 
differential  equation  (1441 


(Y+  1)4>  g2  -  2g  -  g2  +  1  ■  0 
ox’x  ’x  *y 


(A. 45) 


In  this  equation  the  term  <t>QX  which  is  determined  by  the  basic 

field  appears;  it  is  expressed  by  Eqs,  (154).  through  (159)  . 

We  add  the  following  expression. 

It  follows  from  Eqs.  (155)  and  (156) 


x  *My 


V5  r 


Hence,  with  Eq.  (158) 


-x  =  U(a/y2)"2/5(l-(a/2) 


(A. 46) 


and 


(cr/y2)  -M  5/2(-x)’5/2(l-(cr/2))5/2  (A. 47) 

The  following  formulae  will  be  needed,  because  we  shall  use  cr  and 
y  (instead  of  x  and  y)_  as  independent  variables. 

One  has,  from  Eq.  (158) 

dcr/d?«  (10/3)o7/5(o+  (4/3)  j"1.  (A. 48) 
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Then  with  Eqs.  (159) 

3a(x,y)/3x  *  (dc/dcJ  (3~/3x)  * u  1  (10/3) y~4//5o  7//5(a+  (4/3)) 

(A. 49) 

3a (x,y) /3y  =  (da/d ?  (3c/3y)  =  -  (4/3) y-1a(a-2) (o+  (4/3))  1 
It  follows  from  Eqs.  (159)  and  (A. 48)  that 

d7/dT*  (d7/do)(da/dT)  =  a1/5(a/4)-l) 

Then  from  Eqs.  (154), (155)/  and  (156) 

(y+1)*  =  y2(a/y2)1/5(  (a/4)-l)  .  (A. 50) 

vj  X 

As  a  preliminary  hypothesis  we  introduce,  instead  of  Eq.  (161) 

gc  -  u'1  y6/5  7  „<«> •  <A-51) 


Then 

g  -  y’2  (10/3)  y2/5  a7/5(a+  (4/3)) -1  dyQ/da  (A. 52) 

OX 

gQy  =  y^y1/5  [  (6/5)7Q  “  (4/3)  a(a-2)  (a+  (4/3))  1  dy0/da]  . 

This  expression  is  substituted  into  Eq.  (160) .  One  obtains  with 
Eq.  (A. 50) 

li”2y2/5|a1/5(  (a/4)-l)  (100/9)  a14/5(a+  (4/3) )  "2  (dy0/da)  2 
-(20/3)a7/5(a+  (4/3) ) -1  (dyQ/da) 

-1(6/5) y0“(4/3)a(a-2) (a+  (4/3) ) “1dy0/dal 2  }  »0 
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Next  we  set 


-  -  "3/5  ,  . 

Yc  “  0  Y0(<*) 

d70/da  =*  a“8/5("(3/5)y0  +  o(dYQ/da)). 


(A. 53) 


Then,  except  for  a  common  factor,  one  obtains  integral  powers  of 
a  in  the  last  equation 


u"2y2/V6/5{  (100/9)  O  ( <o/4)  -1)  (a+(4/3) )  "2  (-  (3/5)  Yq+  (dYo/da) )  2 

- (20/3) a (a+  (4/3) ) _1 (- (3/5) Yq  +  a(dYo/da))  (A. 54) 

-  ( (6/5)  Yq-  (4/3)  (a-2)  (ff+  (4/3))"1(-(3/5)Y0  +  (dYQ/dV  )  ]2  }  =  0  . 

The  coefficient  of  Y  in  the  last  term  reduces  to 

o 

(6/5) Yq- (4/3) (a-2) (e+  (4/3) ) _1 (- (3/5) Yq 
=  (a+  (4/3))_12 ayQ. 

One  thus  obtains  for  the  whole  expression  (Eq.  (A. 53)) 

U~2y2/5c"1/5(a+  (4/3) )-2  {(100/9) ((o/4) -1) (-(3/5) YQ+  (dYQ/d a))2 

-(20/3) (0+  (4/3)) (-(3/5) Yq  +  a(dYQ/do)) 

«0[2YO  -  (4/3)  (0-2) (dYo/d0)]2  }  -  0. 

Next  we  collect  within  the  braces  the  quadratic  terms  with 
different  powers  of  dYQ/do 


Quadratic  terms : 

(dYo/d0)2I(lOO/9) (0/4)-l)02  -  0(16/9) (o-2) 2 

-  (dYo/do) 2o[o2  -  4o  -  64/9)  -  (dYo/do)2o(o  +  4/3) ) (0  -  (16/3)) 


Linear  terms 


(dy0/da) o[ (100/9) ((a/4) -1) (-6/5) yq  -(20/3) (o+  (4/3)+  (16/3) yq (a-2) ] 

=  (dYo/da)a[Yo(8/3  +  2a) - (20/3) (a+  (4/3))] 

*  (dYo/da)a(a+  (4/3) ) I2Yo- (20/3) ] . 

The  terms  within  the  braces  which  do  not  contain  dy0/da  are 

(100/9) ( (a/4) -1) (9/25) Yq  -  (20/3) (ff+  (4/3) ) (-3/5) Yq  -4 ay*) 

*  ay 2  -  4Y3  +4(a+  (4/3) )yq  -  4oy2 

*  (a+  (4/3) (-3YQ+4) Yq. 

One  thus  obtains 

y"2y2/5a“1/5(a+  (4/3))_1{a  (a- (16/3)) (dYo/da) 2+a (2Yq- (20/3) (dYo/da) 

+  (-3Yq  +  4)yc}  »  0. 

The  solution  which  is  smooth,  even  at  a  *  16/3  (the  limiting 
characteristic!  is  immediately  obvious 

Y  ®  4/3  (A. 55) 

'o 

Combining  Eq.  (A, 51),  (A. 53),  and  (A. 55),  one  obtains 

gQ-  (4/3) y"1  y6/5  a"3/5  (A. 56) 

A  correction  g^  to  gQ  is  computed  from  Eq.  (166) 

<<Y+  D*0!t90x-1)?lx  -  goygly  +  U/2)  m  0  (A’57) 
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One  has  from  Eqs.  (A. 561  and  (A. 49) 


gox  =  -p'2{8/3)y 2/5a'1/5(a+  (4/3)  T1 
goy  -  u"1(8/3)y1/5a2/5(a+  (4/3))"1 
and,  with  Eq.  (A. 50) 

(y+l)<f>oxgox  -  -  (8/3)  ((o/4)-l)  (o+  (4/3)  T1 
(Y+1)<l>oxgox  -1  =  ("  (5/3)  a  +  (4/3))  (a+  (4/3)) 


(A. 58) 


(A. 59) 


First  we  set 

gx  =  y4/5  yx(a)  (A. 60) 

Then  one  has,  with  Eqs.  (A. 49) 

glx  “  w"1d0/3)o7/5(ff+  (4/3))"1(dy1/da) 
gly  =  y"1/5[(4/5)y1  -  (4/3)  (0-2)  (o+  (4/3) ) “1o (dy^/do)  1  . 
Substituting  this  into  Eq.  (A. 5 7) ,  one  obtains 
(-  (5/3)  a+  (4/3) )  (a+  (4/3) )  “1  y_1  (10/3)  a7/5  (o+  (4/3) )  "1  (d^/da) 
-y"1(8/3)y1/5a2/5(o+  ( 4/3) ) -1  y“1/5  [  (4/5) ^-(4/3)  (a-2) 

(a+(4/3))71a(dy1/do)J  +  (1/2)  =  0. 

This  equation  simplifies  : 

(o+  (4/3) )  ”1  u"1  o2/S  |-2  (dYx/da  -  <32/15)7^  +  (1/2)  =0  (A. 61) 
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This  shows  that  with  the  hypothesis  (Eq.(A.57)),  one  obtains  an 
ordinary  differential  equation  in  which  y  does  not  appear. 

Next,  one  introduces 


=  U0  2/^y^.  (A. 62) 

Substitution  into  Eq.  (A. 61)  gives 

|(4/5)Yl  -  2a(dYl/do)  -  ( 32/15) Yl |  =  -(1/2) (c+  4/3)) 
and  finally. 


o(dYl/da)  +  (2/3)Yl  =  (1/4) (o+  (4/3)). 

The  solution  which  is  smooth  at  a  s  0  is  given  by 

Yl  =  (1/2)  +  (3/20)o  (A. 63) 

Thus, 

gx  -yy4/5  o"2/5((l/2)  +  (3/20)a)  (A. 64) 

From  Eqs.  (A. 56)  and  (A. 64) ,  the  following  approximation  is 
obtained 


g  »  y"1y6/5a"3/5[(4/3)+A((l/2)  +  (3/20) o)]  (A. 65) 


with 


A  -  uJ(c/y2>1/5 


f 


In  the  vicinity  of  the  negative  x-axis,  one  uses  Eq.  (A. 47) 
to  evaluate  (o/y2) 
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Next,  formulae  for  the  characteristic  directions  are  developed 
on  the  basis  of  Eq.  (147) .  Parts  of  the  numerator  and  denominator 
of  the  right  side  are  found  in  Eqs.  (A. 58)  and  (A. 59) .  For  the 
missing  contributions  of  g^,  one  finds  from  Eq.  (A. 64)  and  (A. 49) 


y y_1/5o“2/5  { (4/5) ((1/2)  +  (3/20) a) 
-(4/3)  (a-2) (a  +  (4/3) )-1  [-(2/5)  (1/2) 


(A. 66) 

+  (3/5) (3/20) a]| 


This  simplifies  to 


gly  =  u  (16/15)  y~1/5  3/5(a  +  (4/3)  T1  (A. 67) 

Moreover,  from  Eqs.  (A. 64)  and  (A. 49) 

glx  y4/5a"7/5(-(V5)  +  (9/100)a)y“1(10/3)y"4/5a7/5(a+(4/3)  )~1 
glx  =  (_(2/3)  +  0/10)  a)  (o+  (4/3))-1 
Then  with  Eq.  (A. 50) 

(A. 68) 

(Y+1)*oxglx  *  y2(a/y2)1/5((o/4)-1)  (-(2/3)  +  (3/10)0)  (0+  (4/3) ) -1 

Combining  Eqs.  (A. 67)  and  (A. 68)  with  (A. 59)  and  (A. 58)  and  using 
the  definition  of  A,  Eq.  (A. 65) ,  one  obtains 

1  -<Y+  U*ox9x  -  (ct+  (4/3) ) _1 [ (5/3) 0- (4/3) )  +  A( (0/4) -1) ( (2/3- (3/10) ) 0) 

gy  *  y"1  y1/5o2/5(a+  (4/3) ) “1 (8/3) [1  +  (2/5) A] 

Then,  from  Eq.  (147) 

Dx/Dy  *  u  (a/y2)’2/5y“1[(5/8)'a  -  (1/2) 

+  A((c/4)-l)  Cl/4)-(9/80)c] [l+(2/5)A)"1  CA*te3) 
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To  obtain  Do/Dy ,  we  proceed  as  follows.  It  follows  from  Eg.  (A. 46) , 
that  along  a  curve  o*  o (y) 


dx/dy=uy“1(a/y2)”2/5C(4/5) ( (a/2) -1) +( 3/10) (a+ ( 4/3) ) d(loga) /d (logy) ] 

One  obtains,  by  substituting  into  this  equation  the  value  of  Dx/Dy 
for  the  characteristics,  Eq.  (A. 69) 

(5/8)gz.(l/2)^A((?/4)  -p((l/4)-(9/80)a),(4/5)  {(  /2)_1)+(3/10)  (0+(4/3)  )D(10ga)/D(logy) 
1  +  (2/5)  A 


Hence 

(3/10)  (a+  (4/3)  )D(loga)/D(logy)  - 

»  (l+(2/5)A)’1  {(5/8)a~|+A(  (  a/4) -1)  (( 1/4 )-( 9/80  )a  )-<!+(  2/5)  A)  (4/5)  ((a/2)-l)} 
*  (l+(2/5)  A)  ”1{  (9/40)  (a  +(4/3) ) -A(9/320)  (a  +  (4/3)  (a-(28/15))} 

Hence , 

D ( loga )  /D ( logy)  *  (l+(2/5)  A)  "1  (.3/4)  U-(l/8)  A(a-  (28/15) )  ] 

One  may  develop  the  denominator  with  respect  to  A.  Then  one  obtains 

D  ( loga) /D(  logy)  -  3/4  [l-(l/8)  A(a+  (4/3))]. 

Usually,  it  is  convenient  to  consider  a  as  independent  variable. 

Again,  developing  with  respect  to  A,  one  obtains 

D ( logy) /D( logo)  -  (4/3  [l-(l/8) A(a+  (4/3))]  (A. 70) 

Next,  we  derive  an  expression  also  for  D(log-x)/Dloga.  One  finds  from 
Eq.  (A. 46) 

Dlog(-x)  -  (4/5)  d( logy) -(3/10)  (a +  (4/3))  (1  -  (a)/2) )  _1d(loga)  . 


Ill 


/-  •  %<  . 


I- 


Equation  (A. 70)  is  used  to  eliminate  d(logy) .  One  then  obtains 

(A. 71) 

DTlog(-x) ) /D(loga) =  (2/3) [(l-(5/4)o) (l-(a/2) ) _1+(l/5) A(o+  (4/3) ] . 

No  development  with  respect  to  A  has  been  carried  out  in  deriving 
Eq.  (A. 71)  from  Eq.  (A. 70) .  Although  both  equations  neglect  terms 
of  higher  order  in  A,  they  will  give  the  same  directions  in  the  x,y 
plane.  The  same  holds  for  the  quotient  of  Eqs.  (A. 71)  and  (A. 70) 

D(log(-x) )/D(logy)  =  (A. 72) 

(1/2) [ (1-5/ 4) 0) (l-(a/2) ) _1+(l/5) A(a+  (4/3)) ] [l+(l/8) A(a+(4/3) ) ] ~1 

For  the  vicinity  of  the  y  axis,  one  writes 

D(-x)/Dlogy) 

=  (1/2)  [  (1- (5/4) a)  (-x) (I-ct/2) ) _1-x (1/5) A(a+  (4/3)) ] [l+(l/8) A(a+(4/3) ) ]-1 

and  uses  Eq.  (A. 46)  to  express  (-x) (l-(0/2))-1,  which  for  x  =  0 
gives  (0/0) . 

An  analytical  expression  for  the  characteristics  is  found 
from  Eq.  (A. 70) .  In  the  lowest  approximation  in  A,  one  obtains 

y  =  c  c4/3  (A.  7  3) 

where  c  is  the  constant  of  integration.  To  take  first  order  terms 
of  A,  we  introduce 


u  =  log  y 
v  =  log  o. 

This  is  substituted  into  Eq.  (A. 70)  and  the  specific  form  of  A  is 
introduced .  One  obtains 

Du/Dv  -  (4/3)  +(l/6)y2exp(-(2/5)u) Iexp((6/5)v)  +  (4/3) exp ( (l/5v) ]=0 . 
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In  the  first  approximation ,  one  obtains  the  above  results 

■  (4/3) v  +  log  c. 
o 


Next,  one  sets 


u  =  u  +  u. 
o  l 

and  neglects  higher  order  terms  in  u1 .  Then 

DOj/Dv-U/6)  u2(exp(-  (8/15)  -  (2/5)  logc)  {exp ( (6/5)  v)  +(4/3)  exp  ( (1/5)  v)  ]  *0 
Du1/Dv*(l/6)ii2c"2/5{exp((2/3)v)  +  (4/3)  exp(- ( (1/3)  v)  ]  =  0 
=v2c~2/5[(l/4)exp((2/3)v)  -  (2/3)  exp  (- (1/3)  v)  ]  . 

Therefore,  accurate  to  the  next  order 

u  -  log  c  +  (4/3) v  +  u 2c~2/5 [ (1/4) exp( (2/3) v) - (2/3) exp( (-1/3) v) ] 
y  -  c  a4/3exp[y2c"2/5(l/4)o2/3  -  (2/3)a"1/3] 

For  c  large  (y  large)  the  argument  of  the  exponential  function  is 
small.  Developing  the  exponential  function,  one  obtains 

y  -  c  a4/3  +  u2c3/5[(l/4)c2-  (2/3)o]  (A. 74) 

For  the  application  of  the  far  field  conditions  one  needs 
the  derivative  of  g  along  the  characteristics.  The  formula  which 
is  correct  to  the  first  order  in  A  is  Eq.  (177)  .  The  lowest 
approximation  is  simply  by  substituting  Eq.  (A. 74)  into  Eq.  (A. 56) 

ge  -  (4/3)y”1c  6/5 a 
gft  -  (4/3)U-1c  9/20y3/4. 
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Hence , 


D  gQ/D a  *  (4/3)  p"1  c6/5(l  +  0(A)) 

.  =  (4/3)p_1  c6/5(l  +  0(y“2/5)) 


(A. 75) 


and 


D,  VDy  =  y-*=  9/20  . 


Also  needed  is  DlogH/Dy.  One  has  from  Eqs.  (153)  and  (149) 


D  (logH)  /Dy  =  (1/2)  [  (Y+l)  (  3  (4>oxgx)/3x)  -  q  ]/g  . 


(A. 76) 


We  evaluate  this  expression  only  to  the  lowest  order  in  A. 
According  to  Eq.  (A. 59) 


-1 


(Y+  D<l>oxgox  -  -(8/3)  ((o/4)-l)  (c+  (4/3))  . 


Hence  with  Eq.  (A. 49) 


(Y+1)3(<J>  g  )/  =-(8/3)  (1/4)  (g+(4/3-)).-y</4)-l)  y  1  (io/3)  y-4/5a7/5  (0+[4/3) ) 

U  A  OA 


(cr+  (4/3)) 


=  -p'1  (320/27)  to  /y2) 2/5  o(a+  (4/3) )  “3 . 


(A.  77) 


According  to  Eq.  (A. 58) 


(A. 78) 


g, 


oy 


■p-1(8/3)y1/5o2/5(a+  (4/3)’1  =p_1  ( 8/3) y (a/y2)  2/5 (a+  (4/3)  )_1 


Then  with  Eq.  (A. 49) 

p’1(8/3)y“4/5o2/5{(a+(4/3) )_1(l/5)  + 


oyy 


t(2/5)a"1(a+  (4/3) ) "1- (a+  (4/3) ) “2 (- (4/3) a (o-2) (a+4/3) 
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Hence 


9oyyssy"1(8/3)  (0/y3)2/5(a+(4/3))"3{(1/5)°2+(8/15)  +  <16/45> 

- (4/3) (0-2)  [(2/5) <0+  (4/3) )-o]J  (A. 79) 

9oyy  =y”1(8/3) {a/y2)2/5(a  +(4/3)”3[a2_(16/9)  +  (16/9)1 
Then  from  Eqs.  (A. 77)  and  (A. 79) 

(Y+l) 9(0QXgQX)/3x-gyy  =  U~1(a/y2) 2/5 (c+ (4/3) ) ”3 ( 8/3) [-(40/9) a  -  a2 
+  ( 16/9) a  -  16/9]  =  u'1(a/y2)2/5(8/3) (o+(4/3))~1 
and  from  Eq.  (A. 76)  together  with  (A. 78) 

D (logH) dlogy  =  -(1/2)  (A. 80) 

Along  a  characteristic  is  the 

H  =  const  y'1^2 

In  the  expression  for  the  far  field  conditions,  Eq.(80),  DlogH/Ds 
nd  Dg/Ds  appear  in  the  coefficient  of  <j>.  For  large  y, 

DlogH/Dy  'v  y“^ 

-1/4 

and  according  to  Eq.  (75) ,  Dg/Dy  %  =  y  . 

The.  contribution  of  DlogH/Dy  is  therefore  negligible.  Equation 
(80)  is  important  in  the  demonstration  of  Section  IX. 
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APPENDIX  VI 

THE  AXISYMMETRIC  UNSTEADY  FLOW  FIELD  AT  MACH  NUMBER  ONE 

DETAILED  COMPUTATIONS 


The  computations  are  identical  with  those  for  the  two- 
dimensional  case,  except  of  course  for  the  specific  form  of  the 
analytical  expressions. 

The  function  g  satisfies  the  same  differential  equation  as 
before 


(Y+D  4>oxgx  -  29x  -  +  1  -  0-  (A.  81) 

The  term  for  <j>ox,  by  which  the  basic  field  enters  is  found  from 
Eqs.  (186)  through  (189).  We  add,  (from  Eq.  (187)) 

4/7- 
x  =  UY  5 


and  with  Eq.  (189) 

-x  *u  (a/y2r2/7(l-2a) 

and 

WY2)  =  U7/2(-x)-7/2U-2  a)  7/2 
One  obtains,  from  Eq.  (189) 

dc/d^  =  (7/10) o9/7(a+  (1/5)  T1 
and,  with  Eqs.  (186) 

3a(x,y)/3x  =y_1  (7/10)  y_4/7  a9/7  (a+  (1/5) )  _1 
3a (x,y) /3  y  *  - (2/5)y_1  a(2a-l) (a+  (1/5))"1 


(A. 82) 


(A. 83) 


(A. 84) 


(A. 85) 
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- 1  r  ii 


-■  as  t 


It  follows  from  Eqs.  (189)  and  (A. 84),  that 


dT/dc  *  (8/3) a3/? (o-l) . 

Then,  from  Eqs.  (186) , (187) ,  and  (188) 

(Y+1)*ox  =^2^/y2)3/1  (S/3)  (a-D  .  (A.  8 

We  make  the  preliminary  hypothesis 

90  ■  u  1y10//7  yQ  ( a )  .  (A. 8 

Then, 

9ox  =  u"2(7/10)y6/7  o9/7(a+  (1/5)  T1  dyQ/do 

goy  =  y‘1y3/7[U0/7)Yo  -(2/5)a(2a-l)  (a+  (1/5)  “-‘‘dy^da] .  (A. 8 

The  function  gQ  will  satisfy  Eq.  (A. 81)  in  the  highest  power  in  y 
One  has  specifically 

-  2  Sox  -  Soy  -  °’ 

One  obtains,  by  substituting  Eqs.  (A. 86)  and  (A. 88) 
u-2y6/7|a3/7 (8/3) (49/100) a18/7 (o+  (1/5) ) “ 2 (dyQ/da) 2 
- (7/5) a9/7(o+ (1/5) )’1(dy0/da) 

- [ (10/7) - (2/5) a (2a-l) (a+  (1/5) ) ‘1dy0/dc] 2 }  =  0 . 

Next  we  set 


-5/7 

y  Y„ 


dYQ/da  *  a  12//7  [-(5/7)  y  +  a(dy  /do)] 


(A. 8 


Then  one  obtains 


U_2y6/7  a"10/7  |(98/75)a(a-l) (a+(l/5))“2(-(5/7)yo+  a(d70/do))2 

-  (7/5)ct  (o+  (1/5) )  -1  (-  (5/7)  Yq  +  a(dYQ/da))  (A. 90) 

-  [  (10/7) YQ  -  (2/5)  (2a  -1)  (a+ (1/5) ) _1 (- (5/7) Yq+  a (dYQ/da) ) ] 2}=0 . 

The  coefficient  of  Y  in  the  last  term  reduces  to 

o 

(10/7) Yq- (2/5) (2a-l) (a+  (1/5) ) _1 (-5/7) Yq 
=  a (a+  (l/5))‘12Yo. 

One  then  obtains  for  the  whole  expression,  Eq.  (A. 90) 

U~2y6/7  a“3/7(a+  (1/5)  )"2  {  (98/75)  (o-l)  (- (5/7)  Yq  +  a(dyo/dcj))2 
+  (a+  (1/5)) (yq  -  (7/5)o(dYQ/da)) 

-a (2y0  +  (2/5) (l-2o) (dY0/da)]2}  »  0. 

Next  we  collect  within  the  braces  the  terms  with  different  powers 
of  dYQ/da.  Quadratic  terms: 

(dY0/da) 2 [ (98/75) (a-1) a2- (4/25) a (l-2o) 2] 

=(dY0/da)2  [ (2/3)o2-(2/3)a-(4/25) ]  =  (dYQ/dc) 2a (2/3) (a+ (1/5) ) (a- (6/5) ) 
* 

Linear  terms: 

(dY0/da)[ (98/75) (a-1) (-10/7) aYQ- (7/5) (a+  (1/5) a- (/8/5) a (l-2a) yq1 


”(dYrt/da)a(a+  (1/5) ( (4/3) y- (7/5) ] . 


Terms  without  dYQ/do 

(98/75) (o-l)  (25/49) Yq  +  (o+  (1/5) YQ  -  4o/y2 
=  YoI-(10/3)yo  +  1  ](a+  (1/5)). 

One  thus  obtains ,  from  the  terms  within  the  braces 
(o+  (1/5)]  | (dYQ/do)2(2/3) (o-(6/5)  +  (dYQ/do) o ( (4/3) yq- ( 7/5) 
+  Yq[-(10/3)yo  +  1)  }  =0 

The  solution  which  is  smooth,  even  at  a  *  (6/5)  is  given  by 


Y  Q  =  3/10.  (A. 91) 

Combining  Eqs.  (A. 87),  (A. 89),  and  (A. 91),  one  obtains 

gQ  =  (3/10)u_1  y10/7  o“5/7.  (A. 92) 

A  correction,  g1  to  gQ,  is  obtained  from  Eq.  (A. 81)  by  retaining 
only  first  order  terms  in  g1 


((Y+  1)  g  -l)g.  -  g  g  +  (1/2)  *  0. 

ox3ox  *lx  oy^y 

One  has,  from  Eqs.  (A. 92)  and  (A. 89) 

gox  -  -M_2(3/20)y6/7  o~3/7(o+  (1/5) )_1 

goy  -  U-1 (3/5)y3/7o  2/7 (o+  (1/5) )_1 
and  with  Eq.  (A. 86) 

(Y+  "  "(2/5)  (o-l)  (o+  (1/5)  J"1 

(Y+  l)<>oxgox-l-(-( 7/5)o+ (1/5))  (o+d/5))"1 


(A. 93) 


(A. 94) 


(A. 95) 


we  set 


=  y4'/7Y1(a)  (A. 96) 

Then,  with  Eqs.  (A. 49) 

glx  =n"1(7/10)a9/7(a+  (1/5) ) "1dY1/da 

gly  =  y‘3/7[(4/7)71(-2/5)  (2o-l)  (o+  (l/5))“1adY1/da] 

Substituting  this  into  Eq.  (A. 93) ,  one  obtains 
(-(7/5)o  +(1/5)) (o+  (1/5) )_1  y _1 ( 7/10) o9/7 (o+  (1/5) ) "1dY1/da 
-U~1(3/5)y3/7a2/7(a+(i/5))“1[(4/7)Y1  -(2/5) (2a-l) (a+  (1/5) ) ■1adY1/da] 
+  (1/2)  =  0. 

This  equation  simplifies  to 

(A. 97) 

w"1(l/5)(o+  (1/5))"1  o2/7[-(5/2)a(dY1/do)  -  (12/7) ^J  +  d/2)  =0. 

Next  one  introduces 

-  -2/7 

Y1  =  MO  Y]_ 

Substitution  into  Eq.  (A. 97)  gives 

(1/5) [(5/7)Yl  -(5/2)a(dY1/da)-(12/7)Y1]]  +  (1/2) (a+  (l/5)=0 
and  finally 

o(dYl/do)  +  (2/5)  «  (o+  (1/5)). 

The  solution  which  is  smooth  at  a  -  0  is  given  by 
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(1/2)  +  (5/7) a 


(A. 98) 


Thus 

gx  -  uy4/7  a’2/7((.l/2)  +  (5/7) a)  (A.99) 

From  Eqs.  (A. 92)  and  (A.99) ,  one  thus  obtains  the  following 
approximation 


g  -y-1  y10/7  a‘5/7[ (3/10)+A( (l/2)+(5/7)o) ] 


with 


A  -  y2(a/y2)3/7  (A. 100) 

In  the  vicinity  of  the  negative  x  axis  one  uses  Eq.  (A. 83)  to 
evaluate  (a/y2) . 

Next,  formula  for  the  evaluation  of  the  characteristic 
directions  are  developed.  Since  the  differential  equation  for  g, 

Eq.  (144), is  the  same  for  the  plane  and  the  axisymmetric  problem, 
one  obtains  the  same  formula  for  the  characteristic  directions,  viz 
Eq.  (147) .  Parts  of  the  numerator  and  the  denominator  are  found 
in  Eqs.  (A. 94)  and  (A. 95) .  For  the  missing  contributions  of  g^, 
one  finds  from  Eqs.  (A.99)  and  (A. 85) 

gly  *  VY~2/1o~2/1  {(4/7) (1/2) +(5/7) a)  +  (- (1/7) + (25/49) c) 

-  (2/5) (2o-l) (at  (1/5) ) ”1 [ (1/2) (-2/7)  +  (5/7)2al. 

This  simplifies  to 

gly  -  y(24/35)y‘3/7  a5/7(a+  (1/5) )"1.  (A. 101) 


121 


Moreover,  from  Eqs.  (A. 99)  and  (A. 85) 

glx  “ **  y4/?  o"9/7(-(1/7)  +  (25/49)a)y"1(7/10)y'4/7a9/7(a+  (1/5)) 
glx  -  (-(1/10)  +  (5/14) a (a  +  (1/5)  T1. 

Then  with  Eq.  (A. 86) 


(A. 102) 

(r+1)Vlx=lj2(c/y2)3/7(8/3)  (a‘1}  (“(1/10)  +  (5/14)o)  (o+  (1/5)  r1 


-1 


Combining  Eqs.  (A. 101)  and  (A. 102)  with  (A. 94)  and  (A. 95)  and 
using  the  definitions  for  A,  Eq.  (A. 100),  one  obtains 

l-(Y+lHoxgx  -  (cr+  (1/5) )  -1  [  (7/5)  a- (1/5)  +A(a-1)  ( (4/15) -(20/21)  a)  ] 

gy  -  y“1(3/5)y3/7a2/7(+(l/5))“1(l  +  (8/7)A). 

Then,  from  Eq.  (147) 


(A. 103) 

Dx/Dy=  y(a/y2)  _2/7y~1  [  (7/3)  a-  (1/3)  +A(a-1)  (  (4/9)  -  (100/63)  a)  [1+  (  8/7)  A]  _1 


To  obtain  Do/Dy  we  proceed  as  follows.  It  follows  from  Eq.  (A. 82) 
that  along  a  curve  a-  o(y) 

dx/dy  *uy_1(a/y2)"2/7[(4/7) (2a-l)+(10/7) (a+(l/5) )d(logo)/d(logy) ) 

One  obtains,  by  substituting  into  this  equation  the  expression  (A. 103) 
for  the  value  of  Dx/Dy  along  a  characteristic 

IZ/llg."  (I/3.*  th (.?“!)  (Ji/9l-J.l°Q/l3i-?)  a«  (4/7)  (2a-l)  +  (10/7)  (a+(l/5) ) 

(1+ (8/7) A)  D( logo) /D (logy) 


Hence, 
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(10/7)  (o+ (1/5) D(loga)/D (logy) 

=»  (l+(8/7)A)"1{c7/3)o-(l/3)+A(a-l)  (  ( 4/9  H 100/6  3)  o)  -  (l+(8/7)  A)  (4/7)  (2o-l) 
-  (1+ (8/7) A) _1 | (25/21) (a+  (1/5)  -A(100/63) (a+(l/5)  - (23/35) ) | 

Hence , 

D (logo) /D( logy)*  (1  +(8/7) A) ”1 (5/6) [1- (4/3) A(a- (23/35) ) ] . 

Developing  the  denominator  with  respect  to  A,  one  obtains 

D (logo) /D (logy)  *  (5/6) [l-(4/3)A(o+  (1/5))] 

Usually  it  is  convenient  to  consider  a  as  the  independent  variable. 
Again,  developing  the  denominator  one  obtains 

D ( logy )/D( logo)  *  (6/5) [1  +  (4/3)A(o+  (1/5))].  (A. 104) 

We  derive  an  expression  also  for  D(log(-x) )/Da.  One  finds  from 
Eq.  (A. 82) 

d(log(-x))  -  (4/7) d(logy) + (10/7) (a+(l/5)) (2a-l) -1dloga 
Equation  (A. 104)  is  used  to  eliminate  d(logy) .  Then 

(A. 105) 

D (log) -x) ) /Dloga) *  (2/5) [(l-7a) (l-2a) “1+(16/7) A(a+(l/5) ) ] . 

Here  no  development  with  respect  to  A  has  been  carried  out.  The 
quotient  of  Eqs.  (A. 105)  and  (A. 104)  gives 

D(log(-x) ) /D (logy) ■ (1/3) [l-7a) (l-2o) ’1+(16/7) A(a+(l/5) ) ] 

[1  +(4/3) A (a+ (1/5) ) ]-1. 


For  the  vicinity  of  the  negative  y  axis  one  writes 
D(-x)/D(logy)  -  (1/3)  [1-70  (-x)  (1-20  "1+  (16/7)  (-x)A(o+  (1/5))] 

[1+ (4/3) A(a+ (1/5) ) ]_1 

and  uses  Eq.  (A. 82)  to  express  (-x)  (l-2o) _1 . 

To  find  an  analytical  expression  for  the  characteristics,  we 

set 


u  =  log  y 
v  *  log  a 


Then  from  Eq.  (A. 104)  and  (A. 100) 

Du/Dv  =  (6/5)+  (8/5)y2exp(-(6/7)u[exp(  (10/7)  v)  +  (1/5)  exp ( (3/7)  v)  ] 
One  obtains  in  the  first  approximation 


Next  one  sets 


uQ  ■  (6/5) v  +  log  c. 


u  =  u  +  u. 
o  l 

and  neglects  higher  order  terms  in  Then, 

Du-j/Dv-  (8/5)  VI 2  (exp (-  ( 36/35)  v  -  ( 6/7)  logo)  [exp ( (10/7)  v+  (1/5)  exp ( ( 3/7)  v)  ] 
DUj/Dv  -  (8/5) y2c"6/7 [exp (12/5) v)+ (1/5) exp (-( 3/5) v] 

-y2c'6/7[4  exp ( (2/5) v)  -  (8/15)  exp (- (3/5) v) ] . 

Therefore,  correct  to  the  next  order, 
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u  *  logo  +  (6/5) v  +  y2c“6/7[4  exp ( (2/5) v) - (8/15) exp (-3/5) v) 
y  »  ca6/5exp{y2c“6/7[402/5  -  (8/15) a“3/5] } 

For  c  large  (y  large)  the  exponential  function  is  developed.  One 
obtains 

y  -  c  a6/5  +  u2c1/7o3/5(4o  -(8/15)).  (A. 106) 

We  compute  in  the  lowest  approximation  the  derivative  of  g 
along  a  characteristic.  This  quantity  occurs  in  the  application 
of  the  far  field  conditions.  One  obtains  by  substituting  Eq.  (A. 106) 
into  Eq.  (A. 92) 

gQ  =  ( 3/10 )  y"*1  c10/7o 
or 

,o(3A0)»-l«  25/42  YS/6. 

Hence 

DgQ/D0=  (3/10)y'1c1°/7(l  +  0(A)) 

=  (3/10)y"1c10/7(l  +  0(y"1/2))  (A. 107) 

and 

Dgo/Dy«  (l/4)u“1c25/42y-1/6(l+0(y"1/2)) . 

Also  needed  in  the  far  field  condition  Cat  least  until  it 
has  been  shown  that  its  order  can  be  disregarded)  is  D(.logH)/Dy. 

The  characteristic  directions  are  given  by  the  same  general 
formulae  in  the  axisymmetric  and  in  the  plane  flows.  One  obtains 
from  Eqs.  (.185)  and  (149) 


(A. 108) 


D(logH)/Dy  =  (1/2)  [  (y+1)  3  (<J>oxgx)/3x)-y“’13  (ygy)/ay]gy. 

We  evaluate  this  expression  only  to  the  lowest  order  in  A. 

According  to  Eq.  (A. 9 5) 

(Y+  l)<t>oxgox  =  -(2/5)  (o-l)  (o+  (1/5) ) -1 . 

Hence  with  Eq.  (A. 85) 

(Y+l)3(<f>_vgox)/3x  =  -(2/5)  .  u*'1  ( 7/10 ) y“4/7o9/7  (o+ ( 1/5) 

ox  DX  (0+  (1/5) 1 

=  -u”1 (42/125) (0/y2)2/70  (c+(l/5)"3.  (A. 109) 

According  to  Eq.  (A. 94) 

y  goy  =  u”1 (3/5) y10/7o2/7 (o+  (1/5) ) _1  (A. 110) 

Then  with  Eq.  (A. 85) 

y_13(ygoy)/3y  *  u"1(3/5)y“4/7a2/7{(io/7) (o+  (1/5))*1  + 

+  [  (2/7)  a”3,  (o+  (1/5) )  ”1-  (o+  (1/5))-2]  [-(2/5)0(20-1)  (o+d/5)"1] 
y-13(y  goy)/3y*u'1(3/5) (o/y2)2/7(o+  (i/5) ) "3  { (io/7) (o+  (i/5) ) 2 
+  [ (2/7) (o+  (l/5))-a] [-(2/5) (20-1)  } 
y_13  (y  goy)  /Sy-iT1  (3/5)  (o/y 2)  2/7  (o+(l/5) )  “3  (1/7)  {  (10o2+4o+ (2/5) 

+  (2/5) (l-2a) (-5 a+  (2/5))  } 

y_13(y  goy)/3y  -u'^e/s)  (a/y2)2/7(o+  d/5)  )-3(o2+(3/25)o  +  d/25)) 

(A. Ill) 

Combining  Eqs.  (A. 109)  and  (A. Ill)  one  obtains 
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(y+  1)  3(<J>oXgox)/3x  ‘  y”  °{y9oy)/3y 


-  u“1(a/y2)2/7(a+  (1/5)'3  {-(42/125)a-  (6/5 )o2  - (18/125) a- (6/25) )} 


=-u'1(6/5) (a/y2)2/7(a+  (1/5) )“1. 


Hence,  from  Eqs.  (A. 108)  and  (A. 110) 


D ( logH ) /Dy  =  y-1 


H  =  const  y  ^ . 


(A. 112) 


Comparing  Eqs.  A. 112)  and  (A. 107) ,  one  recognizes  that  in  the 
evaluation  of  the  far  field  conditions  D(logHl/Dy  is  negligible  in 
comparison  to  Dg/Dy.  Equation  (A. 112)  is,  however,  important  in 
the  demonstration  of  Section  IX. 


Figure  1.  System  of  Wave  Fronts  for  a  Subsonic  Free 
Stream  Mach.  Number  in  Linearized  Theory. 


Figure  2.  System  of  Wave  Fronts  for  a  Free  Stream 
Mach  Number  One  in  Linearized  Theory. 


Figure  3.  Normal  to  a  Wave  Front  and  Characteristic 
for  a  Wave  Spreading  out  from  the  Origin. 


y-y<> 


Figure  4.  Coordinate  System  for  the  Evaluation  of 
Derivatives  at  the  Edge  of  the  Computed 
Region. 
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Figure  5.  Computed  Characteristics  (Solid  Lines)  and  Their  Analytic 
One-Term  Approximations  (Dotted  Lines)  for  Plane  Flow  with 
the  Free  Stream  Mach  Number  One.  (The  curves  have  been 
made  to  agree  at  their  outer  points) . 


Figure  6.  Computed  Characteristics  (Solid  Lines)  and  Their  Analytic 
Two-Term  Approximations  (Dotted  Lines)  for  Plane  Flow  wit 
a  Free  Stream  Mach  Number  One. 


Figure  7.  One-  and  Two-Term  Approximations  for  Wave  Fronts 
(g  =  20,  30,  40)  in  Plane  Flow  with  a  Free  Stream 
Mach  Number  One, 


Figure  9a.  Same  as  Figure  8,  but  with  a  Different  Starting 
Point  of  the  Computed  Characteristics. 


IS.# 


Figure  9b.  Same  as  Figure  8,  but  with  a  Different  Starting 
Point  of  the  Computed  Characteristics. 


for  Axisymmetric  Flow. 


as.# 


for  Axi symmetric  Flow. 


Figure  13.  Counterpart  of  Figure  8  for  Axisynunetric  Flow 


X 


g(x,y)«g(x0,y0) 


Figure  14.  Coordinate  System  Used  to  Determine 
the  Radius  of  Curvature  of  a  Wave 
Front. 


Figure  15.  Curve  with  Positive  Radius  of  Curvature. 


Figure  16.  Orientation  of  the  System  with. 
Respect  to  the  x,y  System. 


